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ABSTRACT 

The region in north-central Iowa referred to as the Mans on 
disturbed area was investigated with the seismic refraction method 
and the bedrock configuration mapped. The area is approximately 30 
km in diameter and is not detectable from the surface topography; 
however, water wells that penetrate the bedrock indicate that the 
bedrock is composed of disturbed Cretaceous sediments with a central 
region approximately 6 km in diameter composed of Precambrian cry- 
stalline rock. 

Seismic velocity differences between the overlying glacial till 
and the Cretaceous sediments were so small that a statistical program 
was developed to analyze the data. The program developed utilizes 
existing 2 segment regression analyses and extends the method to fit 
3 or more regression lines to seismic data. 

The seismic data was collected on two-mile centers as reversed 
profiles. The data was analyzed and depths and velocities were determined 
and plotted as two separate maps. The resulting bedrock topography includes 
closed depressions, but the circular pattern observed in the subsurface 
data is not reflected in the bedrock topography. Bedrock seismic velocity 
patterns appear to correlate with the Bouguer gravity map developed in 
an earlier study. The Precambrian crystalline rock displayed velocities 
ranging between 3000 m/sec to 5000 m/sec whereas the disturbed sediments 
ranged between 2000 m/sec to 3000 m/sec. 

The study was not able to provide the necessary data to establish 
the origin of the crater. The bedrock configuration map and the velocity 
distribution map do point to another crystalline subcrop west and south 
of the present outcrop. 
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INTRODUCTION 

The identification of granitic well cuttings by the Iowa Geological 
Survey from water well samples taken near the town of Manson, Iowa, 
suggested that an unusual geologic feature was buried at relatively shallow 
depths beneath glacial drift that blankets most of Iowa. A number of water 
wells drilled into bedrock near the wells from which the granitic cuttings 
were taken, revealed that the normal sequence of bedrock formations did not 
exist in the region. Moreover, the cuttings from these wells were highly 
deformed. The Iowa Geological Survey established a boundary enclosing all 
of these wells and labeled the region the Manson disturbed area. 

In an endeavor to explain the occurrence of the granitic rock and 
anomalously soft well-water in this region, the Iowa Geological Survey and 
U.S. Geological Survey conducted a drilling program near Manson in 1952. 

From a bore hole, drilled approximately 6 km north of Manson, 118 m of 
crystalline core was recovered below 28 m (93 ft.) of glacial drift. 

Petrologic analysis showed that it consisted of highly deformed gneiss, 

gneissoid .granite,, granite, and diabase, all Precambrian in age (Hoppin and 

Dryden, 1958). The dimensions of the major components of each rock type 

were also varied. Breccias * with rock fragments as long as 45 cm grading 

down to chloritized microbreccia, made up a large part of the core. 

Weathering alterations occurred as deep as 60 m into the crystalline material 

(Dryden, 1955). Lidiak et al_. (1966) have dated the crystalline rock by the 

g 

Rb-Sr method to be between 1.13 and 1.44 x 10 years. Bunch (1968) found 
that the minerals in the crystalline mass had numerous unusual microstruc- 
tures. Quartz, had well-defined multiple kink bands, and deformation bands 


2 


in plagioclase were observed. 

After an examination of available well logs, both Dryden (1955) and 

Hoppin and Dryden (1958) suggested that the crystalline rock surface is 

2 

nearly flat- topped and covers an area approximately 4 km directly beneath 
the glacial drift. Extending beyond this area, the surface of the crys- 
talline rock falls away in all directions with steeper slopes to the south- 
east. They concluded that the Manson disturbed area has many similarities 
with the characteristics of cryptovolcanic features as outlined by Bucher 
(1936). Bunch (1968) concluded that the numerous structural features in the 
minerals indicated that the Manson disturbed area had undergone violent 
disruption. Irrefutable evidence that the area had been subjected to shock 
pressures that only a meteorite could produce, such as shock-formed glass 
and/or high pressure silica polymorphs, had not been discovered in the thin 
sections . 

A compilation of elevations of the basement rock of Iowa by Yoho (1967) 
shows that at two locations in Ft. Dodge the top of the basement rock is 315 
and 399 m below mean sea level and is typical of north central Iowa. The 
elevation at which the crystalline rock was found north of Manson by Hoppin 
and Dryden (1958) was 351 m above mean sea level. Therefore, the relief of 
the crystalline rock in the Manson disturbed area must be near 700 m. 

A ground water resource study of Webster County by Hale (1955) suggested 
that the regional structure in the Manson area was abruptly broken by 
faulting, forming a roughly circular basin. Figure 1 is a geologic map of 
the Manson area as presented by the Iowa Geological Survey (1969). The 
Precambrian subcrop is shown surrounded by highly deformed lower Cretaceous 
shales. Mississippian limestones and dolomites border the disturbed area on 


Colorado and Dakota Groups. Sandstone and shale with minor limestone. 
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Figure 1. General geologic map of the Manson disturbed area with seismic refraction and borehole 
locations 




GEOLOGIC BOUNDARY ^ " N BOREHOLE • SEISMIC REFRACTION SPREAD 
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the east, and Cretaceous shales and sandstones border on the north, west 
and south. 

In 1953 an aeromagnetic survey was flown over the disturbed area 
(Henderson et al_. , 1963), and the resultant map was incorporated into a 
regional aeromagnetic map of central and southeast Iowa (Henderson and 
Vargo, 1965). Figure 2 is a reproduction of a portion of the 1953 survey. 

A magnetic high with a peak value of 4470 gammas roughly coincides with the 
outline of the crystalline rock as established from borehole information. 
Also prominent on the aeromagnetic map are several other anomalies. An 
elongate anomaly, trending NW-SE and having an intensity of 4605 gammas, is 
centered at a point 6 km northwest of Manson. Smaller anomalies (4115 and 
3620 gammas) are located further from the center of the structure. 

Holtzman (1970) recently completed a gravity survey of the disturbed 
region. He found that the regional effect of the mid-continent gravity 
high, which lies to the southeast of the Manson disturbed area, could be 
minimized with a trend surface technique that was successfully applied by 
Coons, Woollard and Hershey (1967) to similar gravity data. The Bouguer 
gravity map is shown on Figure 3. Holtzman found that a negative Bouguer 
anomaly in the interior of the disturbed area reflected, in part, low 
density brecciated sedimentary rock. Residual maps, calculated from as high 
as sixth degree trend surfaces, suggested to Holtzman that a peripheral 
graben bounds the eastern and southern edge of the Manson disturbed area, 
and that an inner arcuate fault paralleling the graben might also be located 
roughly equidistant between the center of the structure and the graben. 
Holtzman concluded that the gravity data from Manson had an appearance 
similar to gravity data from impact features. With this information at 
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hand the seismic refraction study was initiated to investigate the nature 
and configuration of the upper portion of the Manson disturbed area. 



9 


COLLECTION OF THE SEISMIC DATA 

A total of 107 seismic locations were occupied in the course of the 
investigation. The data were collected at sites loosely fitting a 3.2 x 3.2 
km (2 x 2 mi.) grid to take advantage of the spacing of the county roads in 
the area. In many cases the seismic locations could not be exactly located 
at the grid points because it was necessary to avoid primary roads, farm 
drainage tiles, and electric hi-lines. For identification purposes, the 
seismic sites were numbered along profile lines east or west of a north- 
south trending center line. The center line was drawn 1.6 km (1 mi.) east 
of the west edge of the townships in range 31W, and 11 profile lines per- 
pendicular to the center line were drawn with a spacing every 3.2 km (2 mi.). 
See Figure 1. An example of the identification of a station would be the 
seismic site marked 3-2E. It would be the second seismic site east of the 
center line on line 3. Where a seismic site lies between two profile lines, 
the number is marked on the figure. All grid locations were occupied except 
8-4W and 10-4E. 

The seismic data were collected with a truck-mounted, 24-channel, 
analog, seismic recorder. Single 14 Hz geophones were spaced 32 m apart 
along two 12-channel seismic lines. The low-pass filters on the amplifier 
banks were locked out so that the lower effective frequency was the 14 Hz of 
the geophone response. A high-pass filter of 120 Hz was used throughout the 
field operation. Amplifier gains were adjusted to enhance the first arrival 
breaks.. Timing on the seismic records was marked by timing lines spaced 
0.010 seconds apart. Two independent 100 Hz oscillators were simultaneously 
recorded and later checked against each other, giving a precision of the 
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lines better than 0.001 second in five seconds. 

Two basic geophone spread configurations were employed. At 29 loca- 
tions, most of which were along line 6, a split spread geometry was used, 
and a charge was detonated between them. This geometry produced two sets 
of 12 arrival times for each shot. In several locations on the eastern 
side of the study area in the undisturbed region, the two cables were strung 
perpendicular to each other; in this fashion the presence of a marker 
horizon below the bedrock surface might be used to detect the true bedrock 
dip. Also at 13 other sites, a second charge was detonated at one end of 
the seismic line so that true reversed profiling might be accomplished on 12 
traces. At the 78 other sites, the spread geometry consisted of stringing 
the cables along a straight line, and charges were detonated at each end of 
the seismic line. This geometry produced 24 arrival times for each shot. 

Five to six pound (about 25 NT) charges of 40% strength stick dynamite 
were buried at a depth of approximately 3 m in the glacial drift. The 
charge was placed beyond the farthest geophone at a distance in line with 
the seismic line of from 20 to 60 m and at a distance perpendicular to the 
line of from 3 to 15 m. 

From the seismic recordings the first arrival times and time breaks 
were picked with a 3X power binocular microscope mounted with an eyepiece 
reticule. Three trained observers analyzed several records, and a compari- 
son showed that the arrival times could be picked with a precision of 
± 0.002 seconds. Because the mechanical alignment of the galvanometer and 
timing line optics in the recording oscillograph cannot be accurately 
assured below this level, the overall accuracy of the arrival time deter- 
minations was sufficiently low to warrant the use of one observer for the 
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bulk of the data. 

The raw time values and distances were encoded and digitized. Each 
coding form was sight verified with the original seismic record, and the 
digital records were then mechanically verified with the coding forms. 
Several consistency checks were made on the digitized data. As an example, 
on 24 trace records the time shift and offset distances on traces 1-12 were 
compared to those of traces 13-24, and the difference in ground distance 
between trace 12 and 13 was checked against the standard takeout spacing. 

Several corrections were applied to the time and distance data. A 
relay closure time of 0.001 second was added to each of the raw first 
arrival times. The horizontal distance from shot point to geophone for 
each arrival time was corrected for the distance the shot point was offset 
from the seismic line. The relief along the geophone spread was usually 
less than 5 m; only in one instance was the relief as great at 15 m. 

Because of this low relief, time corrections due to topography were not 
calculated. 

A reconnaissance altimeter survey of many of the gravity stations 
established by Holtzman (1970) confirmed his published ± 1.2 m elevation 
accuracy. For each seismic location an altimeter loop was tied from the 
seismic shot points to the nearest gravity station. The horizontal datum 
plane for each shot, from which the depth calculations were made, was set 
at an elevation, E, such that E = E shot - C/2, where E s ^ ot is the elevation 
of the top of the shot hole and C is the depth to the center of charge. 
Therefore the accuracy of the elevations of the refracting horizons is 
within ± 2.5 m. 
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STATISTICAL METHOD 

In the development of the refraction method, it was immediately 
realized that there was some reproduceable linear relationship between the 
source-to- receiver surface distance and the time of arrival of the seismic 
energy. The exact linear relationship, as inferred from geometrical optics, 
is often expressed as a line segment on a time-distance plot of the data, 
and the least squares method has been used to 'best fit 1 these line segments 
to the data. 

Background 

The statistical implications of the 'best fit' regression analysis have 
not, until recently, been appreciated. Steinhart and Meyer (1961 ) have 

5 

reviewed the use of statistical uncertainty in seismic refraction, and they 
have concluded that the results of seismic refraction investigations have 
seldom been presented with estimates of uncertainty. 

Because the error of measurement of distance is normally several orders 
of magnitude less than the error of the measurement of time, as a reasonable 
approximation, the error associated with the data may be considered to arise 
only from time. The measurement of distance is then taken to be without 
error. From this, all deviation from the time value predicted by this 
linear relationship are thought to arise from the error of the measurement 
of time only. 

The principle of least squares states that the best estimate of the 
true time value at a distance is one that minimizes the sum of the squares 
of the time deviations from their estimates. The wave slowness and time 
intercept, determined with least squares regression, are used as 
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coefficients of the regression line of best fit. But by the admission of 
error, there cannot be considered an 'exact' functional relationship between 
time and distance, and so, the measurement of the coefficients is also 
subject to error. 

No assumption as to the distribution of time deviations at a distance 
need be made. The principle of least squares does not imply the normal 
distribution of errors, only that the estimated parameters are unbiased and 
have minimum variance. Borcherdt and (Healy (1968), in fact, suggest that 
for crustal refraction data, the distribution is not normal but is more 
peaked than normal. Regardless of whether or not it is, the central limit 
theory states that the distribution of the sum of a number of small errors 
will approach a normal distribution when the number of measurements grows 
large. If no assumption as to the distribution of time error at a particular 
distance is made, or no measure of the distribution is attempted, only the 
determination of the slope and time intercept parameters with the method 
of least squares may be made. Ultimately the time deviations at a distance 
are not as important as the ability of the model to fit the situation. 

The method of least squares was first used to choose the line that 
'best fit' the data over a range of distances chosen by the seismologist. 
Steinhart and Meyer (1961) presented the uncertainty of the wave slowness 
and time intercept for a line segment under the assumption that the error 
of time was normally distributed. Statistics may further be used to choose 
the 'best fit' combination of line segments for the entire range of dis- 
tances over which the refraction data was collected. This problem has not 
been attacked in the field of seismology. The problem has only recently 
been investigated by those in the field of statistics. 




14 


Sprent (1961) considered the determination of the critical value of 
the independent variable when there were known to be two connected regression 
lines. Quandt (1958b) considered the estimation of parameters when the 
regression system obeys two separate regimes. He concluded that some 
previous determination must be made as to the number of regressions present 
in the system. Several difficulties were encountered by Quandt (1958a) when 
he attempted to establish statistical tests to determine this number. 

Robinson (1964) investigated the change-over value for polynomical regres- 
sions. Hudson (1966) presented the general solution to the estimation of 
parameters and change-over values for both linear and higher order poly- 
nonomial regressions. 

Two computational schemes, following Hudson's guidelines, have been 
developed to estimate the number of regressions which will fit a set of 
data. Bellmen and Roth (1969) have used a dynamic programming approach, and 
McGee and Carleton (1970) have applied the method of cluster analysis. 

Recently seismologists have approached the determination of join point 
(or critical distance) by using selective data rejection. Iyer ejt a]_. (1969) 
used an iterative scheme which rejects data with high time deviations. Seg- 
ments of the data are selected to be fitted, and at each successive itera- 
tion, data points which deviate from the calculated regression line by some 
specified level are rejected from subsequent calculations. The curve 
fitting method formulated by Kail a and Narain (1970) differs from the method 
used by Iyer et al_. in that weights are assigned to each data point. The 
weights are assigned on the basis of the normal probability function for the 
time error. Then the regression for each line segment is calculated using 
all the data points weighted for that particular regression. Because of slow 
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convergence, however. Kail a and Narain found it necessary to revert to data 
rejection in the vicinity of the critical distance at higher iteration 
levels. 

A method is presented here that makes no supposition as to the quality 
of the data and takes advantage of several characteristics of seismic 
refraction data. Because shallow seismic refraction data may be modeled by 
a series of line segments, many aspects of the general solution of the 
multi segmented regression problem, as presented by Hudson, may be omitted. 
The parameters for each line segment, wave slowness and time intercept, are 
calculated, and under the assumption that the distribution of time values 
at a distance are normally and independently distributed, the standard 
deviations of these parameters are found. An algorithm, utilizing Hudson's 
2 segment regression analysis, is presented here which extends his method 
to fit 3 or more regression lines to seismic refraction data. The statis- 
tical uncertainties of the parameters for each line segment are used to 
estimate the uncertainty in the determination of depth. 

Unconstrained Join Type 

Let the n time and distance pairs form a series in ascending order, 

and let there be only one time value for each distance, as [(t-j, d^) 

(t n ,d n )]. Let the t^ be the dependent variable, and let the d^ be the 
independent variable. Now the set of time-distance pairs may be partitioned 
into m sections, where j = l,...,m, and a simple linear regression equation 
may be applied to each section of data. Let one of these sections be 
composed of a series of s pairs of values t(tg»dg),.. .,(t^,d^)]. Here 
s = h-g+1. The model equation for the jtti_ section of the data is 
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t . = a- + 3 .d. + £.. 
l J j 1 lj 


i = g h 


The line parameters a. and 3- are the time-intercept and wave slowness, 

J J 

respectively. The random variables e.. are normally and independently 

* J 

? 

distributed with a mean of zero and variance a. . Let the sample estimate 

J 

for the time intercepts be E(a.) = a. and the sample estimate for the wave 

J J 

2 2 

slowness by E(3J - b.. The estimator of the variance is E(a. ) = s. . 

J J J J 

The estimator of each e.. is zero, so the m regression equations are 

. « 

determined by calculating (a^, b-j,..., a m> b m ). Statistical weights could 
easily be included here, but for the present discussion they have been 
omitted. 

The line segments should be continuously connected at points D. which 

J 

are positioned between the end elements of the partitions. There must also 
be at least 2 points of data on each segment, then the time estimates are 


t. = a, + b n d- 
i 1 li 


for 0 <d-j<d 2 <D-| 


ti = a 2 + b 2 d. 


for D 1 <d 1 -<d 1 . +1 <D 2 


t- = a . + b - d. 
i j j l 


for 


t • = a + b d. 
i m mi 


for D m <d n -i <d n 
m n-l n 


For the jth section of the time-distance series, the deviation of each time 


value, tj ,from it's estimation, t. ,is 


e. . = t. - a . - b .d. 
ij l j j i 
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Temporarily let the jth^ section of the data be considered separately, 
as [(t^.d^) i = g,...,h]. The residual sum of squares from Equation 3 is 

s p 

RSS . = E e/ (4) 

J k=l K 


where the length of the jjth^ section is s. 
calculated. 


The following sums must also be 


Sd. 


s 

E d, 
k=l k 


St. = E t, 
J k=l K 


Sdd . = £ d, 

J k=l * 


stt. = E t.‘ 
J k=l K 


(5) 


Sdt . 

J 


k=l K K 


Steinhart and Meyer (1961, Eq. 1 and 2) have presented the solution of 


the line parameters, a. and b.. They are 

J J 

Sdd. St. - Sd .Sdt . 
a , = i J i — _ — 1 


j ssdd. - (Sd .) 2 

J J 


( 6 ) 


and 


sSdt. - Sd .St . 

b = J iL -4 

J sSdd . - (Sd-T 

J J 


(7) 


Now the variance of the wave slowness is 


a /(b) = RSS j /t(s-2)SSd j ] 


(8) 


SSd . = E (d - <d>) : 
J k=l K 


where 
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The angle brackets indicate average. Steinhart and Meyer (Eq. 13) have 
shown that the variance of the mean value of the regression line at any 


distance, d,is 

o/(d) 




(d-<b» 2 ~ 
SSd . 


( 9 ) 


At the time intercept, d = 0, the variance is 


Oj 2 ( a) = (RSSj/s-2) (1/s + <d) 2 /SSd j ) 


( 10 ) 


Substituting Equation 8 into Equation 10 and using the identity 
s s 

SSd. = z d 2 - ( Z d, ) 2 /s, then 

J k=l K k=l K 

(ID 

Oj 2 (a) = [o. 2 (b) Sddj ] /s 


Equations 4, 8, and 11 may be applied to any section of the data. 

The join points D^,...,D represented on the model given by Equation 2 
are not constrained to a particular value but may be located anywhere 
between the uppermost distance d^ of one segment and the lowermost distance 
d 1>+ .| of the next segment. Let this type of join, where the data are 
separated, be called s-join. The total residual sum of squares is 

m 

RSS t = Z RSS. (12) 

T j=i J 

Consider a 2-segment regression (r = 2). The computation necessary to 
find the best statistical fit of 2 lines fitted to the data proceeds by 
passing a line through [(t^,d^), (tg.dg)] without error and regressing a 
line through [(t^,dj) ( t n ,d n ) ] . If the intersection of these 2 lines 
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is between d 2 and d^, then these are a valid pair of regression lines 
because the regression lines join where the data have been separated. The 
total regression sum of squares RSSy = RSS-| + RSS 2 for this combination of 
lines. If instead the intersection is not between d 2 and d 3 this combina- 
tion of lines is rejected as a possible 2-segment regression candidate. 

Next the first (3) pairs and the last (n-3) pairs are fitted. Again, if the 
intersection is found to be between d^ and d^ we have an acceptable candi- 
date. The procedure is continued, testing each possible combination of 
data, to the last combination which is [(ti,d^),...(t 2 ,d 2 )] and 

[(t^ n ), (t_,d )]. Of these candidates the best statistical fit is the 

n-1 n- r n n 

candidate. which has the smallest RSS-p For r = 2 the total number of 
combinations that must be tested is (n-2). 

Hudson (1966, p. 1104) has shown that for those candidates that have 
not been rejected, the estimates of the regression parameters are valid. 

After the best statistical fit has been found, the coefficients (a^,a 2 ,b-|, 

p p 

b 2 ) are calculated from Equations 6 and 7 and their variances [o^ ( a) ,o-j (b), 
2 2 

a 2 (a),a 2 (b)] are calculated from Equations 8 and 11. 

If r>2 then the same basic computational plan would follow, and if all 
the join points were s-join, the overall least squares solution would be 
the best fit of r line segments found as local least squares solutions from 
all possible combinations of the data points. 

Constrained Join Type 

Although the data have been searched for all s-join candidates, the set 
of all possible linear regression fits is still not complete. The model 
regression equations in Equation 2 must also be tested for the case where 
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the join points are located on time ordinate fixed at the data point 
distances. This type of join is here called c-join because, for each pair 
of regression lines, one data point is common to the solution of each line. 
The system of model equations, with the constraint that the lines join 
somewhere on ordinate corresponding to d. , has been solved with the use of 
a Lagrange multiplier. 

As before, the details of the procedure will be presented for r = 2. 
The appropriate linear constraint is that the lines join at some D. 


a-j + b-jD = ^ + ^2° 


(13) 


Again let the time-distance pairs be [(tpd^) (t n ,d n )], and let 

the c-join be located at D = d . Before the calculations of the residual 

Y 

sum of squares and the coefficients of the line segments are made, a 
number of sums must be found. 


Y 

Sd, = Z d u 
1 k=l 

• Y 

Sdd r = Z d u 
1 k-1 K 


n 

Sd 9 = Z 
L k=y+l 


Sdd 9 = Z d 
L k=y+l k 


(14) 


Following Hudson (p. 1126) the terms <p and \ p are 


and 


where 


<p = a-|-a 2 + (b 1 -b 2 )D 

Sdd, - 2DSd, + yD 2 Sdd. - 2DSd ? + (n- Y )D 2 
m = + 

r o)-| 0^2 

to-j = YSdd^ - ( Sd-j ) 2 

u>2 = (n-Y) Sdd 2 - ( Sd£ ) 2 


(15) 

(16) 
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* 

Let RSSj be the total residual sum of squares for two line segments 
connected with a c-join, and from before, let RSSj and RSS 2 be the uncon- 
strained residual sum of squares, then 

RSSj* = RSS 1 + RSS 2 + <p 2 /ip (17) 

^ ^ ^ ^ 

Finally, let (a-| , b-j , a 2 , b 2 ) be the coefficients of the line segments 
joined by c-join, and let (a^, b^ , a 2 , b 2 ) be the unconstrained coefficients 
as determined from before. The new coefficients are: 

a-| = a.| - ( ) [Sdd-j - DSd^ ] 

b-,* = b 1 + (4>/^co-, ) (Sd 1 - yD] 

08 ) 

a 2 = a 2 + ^ Sdd 2 " ^Sd 2 ] 

b 2* = b 2 ' 0/^ 2 )[Sd 2 - (n-y)D] 

A line segment with these coefficients does not pass through its point of 
mean time and distance. 

The computational scheme could proceed from this point by testing the 

★ 

data from y = 2,...,n-l and finding the minimum RSSy , as had been done for 
the s-join candidates, but several short cuts may be taken. The complete 
set of combinations with s-joins may be used to reduce the number of 
combinations with c-joins. Hudson (1966, Theorem 3) has investigated the 
case for r = 2. Each of these combinations may be rejected and not 
calculated if they pass one of the following 5 tests: 

Test 1. At each s-join candidate, say at i = £, the c-join candidate 
need not be calculated at y = £ + 1. 
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Test 2. If at any i, RSS-j + RSS 2 > Min RSSj (Minimum RSSj) , even if 

the join point was not between the end members of the 2 sets of data, the 

c-join candidate need not be calculated at y = i, i+1. 

At this point the static tests for the elimination of possible c-join 

* 

candidates have been exhausted, and the calculation of RSSj must be made 

* 

at c-join locations from y = 3 n-2. After each calculation of RSSj 

through, the next 2 dynamic tests should be made to further reduce the 

number of possible c-join locations. 

* * 
Test 3. Let a RSSj be calculated and from these let a Min RSSj be 

* 

found. If at any i, RSS-j + RSS^ > Min RSSj , even if the join was not 

located between the end members, then the c-join candidate need not be 

calculated for y = i, i+1. 

* 

Test 4. Let a RSSj be calculated at some y and that the mean dis- 
tances for the left and right line segments are found as (d.|) = Sd^/y and 
(d r ) = Sd 2 /(n-y). If the join point D from the s-join calculations is 
(d-|) < D < (d r ) then RSSj need not be calculated at any y such that 
< d !> <. d y < <d r >. 

Finally a test at y=2 and y = n-1 is presented. 

Test 5. At y = 2 and y = n-1, regression lines may be fitted and 
residual sum of squares, RSS^ and RSS n _.j, may be calculated on the data 

sets [(t 2 ,d 2 ),...,(t n ,d n )] and [ ( t ] ,d 1 ) (t n-1 ,d n _ 1 ) ] , respectively. 

Then a second regression line may be drawn through (t-pd-j) and the inter- 
section of the first regression line with d 2 (without error). Similarly a 

second regression may be drawn through (t n ,d n ) and the intersection of the 

★ 

second regression line with d n _-|. If RSS 2 > Min RSSj, then RSSj at y = 2 
need not be calculated. And if RSS n _-| > Min RSSj, then RSSj at y = n-1 
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need not be calculated. 

At this point for r = 2, all (n-3) calculations have been completed 

for the s-join candidates, and of these, the 2-segment candidate with Min 

RSSy has been retained. There are (n-2) possible candidates with c-join, 

and hopefully, many of these have been rejected by the tests. From those 

* 

that were calculated, the c-join candidate with Min RSSy has also been 

retained. There are then a total of (2n-5) possibilities for r = 2. The 

final least squares solution is either the minimum s-join candidate or the 

minimum c-join candidate, whichever has the smaller residual sum of squares. 

As before when r>2, the same basic computational scheme would be 

.followed for the search for c-join candidates. It is not surprising that 

the total number of possible candidates 'grows rapidly for problems involving 

many segment models or large sets of data, although many may be rejected by 

the tests. When n = 24 and r = 2 the number of possible candidates is 43, 

and when n = 24 and r - 3 the total is 841. For greater n or r the number 

of possible candidates is extremely large. For example, when n = 100 and 

g 

r = 6 the total number of possible candidates is near 1.9 X 10 . 

An Algorithm 

Seismic refraction measurements provide data which are sufficiently 
linear to allow certain liberties in the curve-fitting process. The first 
of these is to suppose that some reasonable estimate may be made of the 
maximum number of permissible line segments present in the investigation. 

The estimate does not need to be known, but that only some upper limiting 
number may be chosen. 

This decision though must be tempered with the inherent tendency for 
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higher-numbered regression systems to fit the data better than lower- 
numbered systems. To be guided in this direction, one must choose a value 
high enough so that no line segment is missed, but not so high that the 
curve is ‘over-fitted 1 and the solution is physically impossible. This 
would be the case when the fitted solution contains velocity reversals. 

The second criterion then is that starting from the origin, the velocity of 
each line segment must monotonical ly increase. In practice, this is a 
valuable guide to the selection of the best least squares solution. 

In a 2 segment linear regression model, the residual sum of squares may 
be plotted against join point distance. In the limit as n °° the residual 
sum of squares form a curve, and the minimum point on this curve corres- 
ponds to the distance at which the two regression lines intersect. For a 
3 segment linear regression as n -► °° the best least squares solution is the 
minimum point on a surface of residual sum of squares. This surface is 
formed on 2 axes which correspond to the D-j and join point values and a 
third axis which is the residual sum of squares. The anology is similar 
for more line segments. If this r-dimension surface is reasonably simple, 
the minimum point may be reached by making a series of traverses across the 
surface, each of which is accomplished by varying the values of only one 
join point. Each traverse is equivalent to the 2 segment regression problem 
considered earlier. 

To evaluate an r>2 regression system, a number of points along the 
data must be chosen to effectively lower the number of possible candidates 
to be searched at any one time to (2n-5). These test points should be 
chosen as near as possible to the suspected join points, but it is prefera- 
ble to choose them slightly below the suspected join point. In this way no 
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data points from the next higher line segment enter the calculation. The 
procedure advances by searching for 2 segment joins in successive sets of 
data. Each successive set of data has an upper bound which has been set by 
a test point and a lower bound which has been set at the join point found 
in a preceding join search. 

As an example, suppose there are 3 suspected line segments to be 
fitted to time-distance pairs where [(t^,d.) i = l,...,n]. There must be 1 
test point chosen at or slightly below the suspected join between the second 
and third line segments, say at i = T. The pairs [ ( t^ ,d^ ) i = 1,...,T] are 
first searched for the join between the first and second line segment. Let 
the solution for this set of data be an s-join candidate between i = s and 
i = c + 1 , as d^cD-j <d ^ + -j - The data [(t^d^) i = ?+l,...,n] is now searched 
for the join between the second and third line segment. Suppose the 
solution here is a c-join candidate such that . The final 

solution is composed of 3 line segments joined at and D 2 , and the total 
residual sum of squares is equal to the sum of the residual sum of squares 
from each line segment. 

If the solution in the first search had been a c-join candidate the 
data in the second search would be scanned from [ ( t^ ,d^ ) i = £,..., n]. Or 
if a join point had not been found in the first search the data in the 
second search would have been scanned from [(t^,d.) i = l,...,n], and only 
a 2 segment solution would have been found. The computation for each search 
would proceed accordingly. 

For an n-segment linear regression a vector of test points must be 
chosen as (T-j ,. . . .T^) and the algorithm proceeds with (n-1) levels of 2 
segment join search. Because the number of possible candidates to be 
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searched at each level of join point selection is restricted to the search 
for a 2 segment fit, the actual computation is greatly reduced. 

Examples 

An example of a set of data from a typical shallow refraction shot is 

presented in Figure 4. Here [■( , d . ) i = 1,...,21]. A linear regression 

with all 21 data points yields a line with a velocity of 2682 m/sec and a 

? 2 

residual sum of squares of 0.47325 x 10 sec . The best s-join candidate 

from a 2 segment regression analysis is superimposed on the data in the 

figure. The join point is located at a distance of 421 m and the total 

-3 2 

residual sum of squares is 0.13222 x 10 sec . All possible c-join 
combinations were eliminated by the 5 tests, so the s-join candidate is the 
best overall 2 segment regression fit. 

The best 3 segment regression is presented in Figure 5. A test point 
at dg was chosen, and a 2 segment analysis of [(t^,d^) i = 1,...,8] revealed 
that the best fit was a c-join candidate joined on the time ordinate 
corresponding to d^. The second level of data then consisted of [ ( t^ , d^ ) 
i = 2,..., 21]. The best statistical fit here was an s-join candidate with 
a join between d-j 2 and d-jg. To be sure the test point was correctly chosen, 
the set [(t.j,d.j) i = 1 , . . . , 1 2 ] was analyzed. This time the best fit was a 
c-join candidate with join point at dy. This third scan of the data 
protected the overall solution from an inaccurate choice of test point. 

Here the test point of dg was too low and should have perhaps been d^* The 
overall 3 segment regression fit has a total residual sum of squares of 
0.44170 x 10 ^ sec^. A borehole adjacent to the shot showed that the 3 
segment fit compared favorably in velocity and depth with glacial drift, 
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shale, and limestone. 

These two examples show that two valid regression models may be fitted 
to the shallow seismic refraction data. At other times only a 2 segment 
regression model may be acceptable. This situation might arise if the 3 
segment regression model contained a velocity reversal between two 
adjoining line segments. In the first case additional information must be 
used to discriminate between the two cases. This information is acquired 
from a knowledge of the geological conditions in the area. 
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ANALYSIS OF NEAR-SURFACE VELOCITIES 

Multi-layered earth models of any dimension, whether they are based 
on bedrock structures measured in meters or crustal structures measured in 
kilometers, may be developed from seismic refraction data only when velocity 
measurement of the uppermost layer has been made. On many scales the earth 
may be approximated by a model bounded by a series of planar surfaces. The 
compressional wave velocity of this uppermost earth material may, in theory, 
be calculated from the reciprocal of the slope of a line segment that 
passes through the origin of the time-distance graph of the data. 

Selection of Data 

Many seismic refraction investigations are hampered by a lack of 
knowledge of the near-surface geologic environment near the shot and under 
each geophone. In shallow refraction investigations the plot of some 
records reveals that perhaps one or two data points nearest the shot fall 
markedly below the first fitted line segment. They may represent random 
scatter about the first line segment, or they may represent travel times 
from direct waves passing through a shallower layer. In either case, the 
first layer is not detected because there is a lack of data from each 
record due to the fixed sampling intervals along the distance axis at which 
time was measured. By combining the data from all records the lack of data 
at near distances on individual records is overcome, and the maximum times 
in each interval may be used to estimate the minimum velocity of the first 
detected layer of a generalized model for the entire study area. 

In the Manson area, line segments drawn through the first several time 
distance points on each shot record did not pass through the origin, so the 
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entire set of first arrival time-distances pairs was plotted. The field of 
points is shown on the scatter diagram (Figure 6). The locus of points 
near the time axis clearly demonstrates that the first layer line segments 
do not pass through the origin. Although the lower border of data scatter 
is not well marked, the upper border is strikingly smooth. 

In an endeavor to better understand the nature of this upper border, 
the distances of all time-distances pairs were partitioned with a class 
interval of 5 m, as 0-5 m, 5-10 m, etc., and the corresponding time values 
in each interval were scanned for the maximum time in each interval. A 
plot of these maximum time values was found to be quite regular even to a 
distance of 500 m from the origin. The maximum time at this distance was 
0.275 sec. The scatter of the maximum times beyond 500 m is probably due 
to the smaller sample sizes in each interval at these farther distances. 

It was found that the locations of these time-distance pairs which had 
a maximum time were uniformly scattered over the entire area. The original 
records, from which these time values were taken, were then reexamined to 
see if these times might have been taken from traces where the refracted 
wave was obsured or where the arrival time was incorrectly measured. 

Because of the large size of the shot charges, first arrival onsets were 
quite strong even at the farthest recording distances. It is unlikely 
that the maximum time values could be attributed to poor record quality or 
to second arrivals from later phases, particularly for distances less than 
500 m. 

Statistical Analysis 

Once that it had been established that the maximum time values were 
valid refraction arrival times, these maximum time-distance pairs were 
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statistically fitted with line segments as if they were data from a single 
seismic record as described previously. By this method a portion of the 
set of all possible combinations of linearly regressed line segments was 
searched for the one fit which had the minimum residual sum of squares. 

The first 100 time-distance pairs were fitted with 1, 2, 3, and 4 line 
segments. The results are tabulated in Table 1. The first line segment 
from the 4 segment model best fits the maximum time values at the nearest 
distances. Models with a greater number of line segments suffered from 
velocity reversals in the first and second line segments, probably due to 
random deviation of the maximum time values and accentuated by the fewer 
number of points represented on the first fitted line segment. 


Table 1. Statistical fits of first 100 maximum times 


Model 

Points 

Velocity 9 

Residual 
sum of. 
squares 0 

1 layer 

100 

1919 

0.24530 

2 layer 

53 

1771 



47 

2138 

0.12464 

3 layer 

12 

1462 



50 

1844 



38 

2203 

0.10938 

4 layer 

7 

1258 - 



20 

1749 



40 

1863 



35 

-2200 

0.10772 


a In m/sec. 
b In 10" 2 sec 2 . 
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From the analysis of the data collected at Manson, it is inferred that 
the minimum velocity of the first detected layer is near 1250 m/sec. A 
line with a slope corresponding to a layer with velocity of 1250 m/sec may 
not be made to pass through the origin and still remain within the data 
scatter. Therefore, the first layer has not been detected, for if it were, 
the line with a velocity of 1250 m/sec would pass through the origin. 

Several investigators (Hobson, 1969; Staub, 1969) have suggested that 
at very short distances, the first detected energy arrives at the geophone 
from a ground-coupled air wave when the velocity of the surficial layer is 
less than the velocity of sound. Although the detection of this air wave 
was not an object of the present study, all depth calculations in this study 
include a first line segment which is attributed to a ground-coupled air 
wave. To avoid confusion the first statistically fitted layer has been 
called the first layer. The depth to the '"air" - first layer' interface 
has been found to range from 1 to 2.5 m. 
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ANALYSIS OF DEEPER VELOCITIES 

The velocity determinations from the statistical line-fitting analysis 
of individual seismic records may be combined to indicate the variability 
of the velocities of various seismic refractors. The utility of the 
shallow seismic refraction technique rests not on the analysis of each 
seismic record, but on how well seismic refractors may be traced below 
many seismic stations, and how well the depths to the seismic refractors 
compare with lithologic changes at depth, as indicated from borehole data. 

Statistical Analysis 

Nineteen seismic spreads were placed as near as possible to boreholes 
where adequate well logs were available, so that a correlation between 
seismic refractor layers and lithologic units could be established. These 
records were plotted and analyzed with the statistical line fitting method 
mentioned previously. Test points were chosen as near as possible to the 
join between suspected second and third line segments. Each seismic record 
was fitted with 2 and 3 line segments. 

The uniqueness of the statistical fit depends on the choice of test 
points. In an attempt to improve the 3 line segment fits, each record was 
analyzed with a variety of test points. The final statistical fit of the 
data for the 3-layer model does not seem to be very sensitive to the loca- 
tion of the test point. Experience on artificially generated 3-layer data 
has shown that if the test point is chosen between the middle of the second 
line segment and the join point between the second and third lines segments, 
a correct solution is found. 

The statistical fitting method was also measured by how well it 
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eliminated the number of calculations of 2 segment regressions. For any 
particular set of n time-distance pairs, the search for the best fit of 2 
line segments joined between data points (s-join) was necessarily applied 
to all (n-3) combination of data points. The robustness of the 5 tests 
used to eliminate the calculation of the regressions of 2 line segments 
connected with a c-join was measured by the percentage of calculations that 
had been eliminated. Of the (n-2) possible combinations to be tested for 
c-join, approximately 90% could be eliminated from consideration for any 
one of the 5 tests. The rejection rate for each of the 5 tests was approxi- 
mately: test 1 - 0%, test 2 - 50%, test 3 - 40%, test 4 - 2%, and test 5 - 

8 %. 

After the validity of the statistical method had been tested, rough 
correlations were made between model velocities and the borehole litholo- 
gies. Test points based on those correlations were then picked on plots of 
seismic profiles without well log control. These seismic records were then 
analyzed for both 2 and 3 line segments. In no instance could 4 line seg- 
ments be satisfactorily fitted to the data because no more than 24 data 
points were available for each statistical analysis. After the best 
statistical fits had been made it was necessary to decide, for each record, 
whether the 2 layer or 3 layer model best fit the geologic conditions. 

Velocity Distributions 

The final results are plotted on Figure 7. The frequency polygons for 
the first layer velocities, second layer velocities, and a combination of 
first, second, and third layer velocities are plotted with the distance axis 
divided by a class interval of 50 m/sec. The results are plotted at the 
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mid-point between the ends of the class intervals. 

The narrow distribution of first layer velocities is attributed to 
glacial drift. Much of the velocity of the glacial drift is between 1800 
and 1850 m/s ec. It is significant to note that the distribution about the 
main peak is symmetric. This indicates that there is only one type of 
material measured in the distribution. The presence of a second material 
with a slightly higher velocity, for example the disturbed shale, would 
skew the distribution to the right. The presence of a second material with 
a slightly lower velocity would accordingly skew the distribution to the 
left. The shape of the main first layer distribution compares favorably 
with the velocity distribution of other glacial drifts as presented by 
Hobson (1969), Lennox and Carlson (1969) and Staub (1969). The several 
first layer velocities above 3600 m/sec were from seismic stations on the 
eastern portion of the area in the undisturbed area where the glacial drift 
was known to be quite thin. 

The frequency of velocities between 2250 and 2300 m/sec on the second 
graph is the highest on the plot. The distribution is not nearly as sharp 
as the glacial drift distribution in the first layer but it represents a 
portion of the disturbed shale velocity distribution. Superimposed on this 
graph are the velocities of the undisturbed bedrock layers which include 
undisturbed limestones and dolomites ranging in velocity from 3000 to 7000 
m/sec and the velocity of the crystalline bedrock which ranges from 3000 to 
5000 m/sec. 

As an example of the irregularity of the bedrock velocity, at station 
6-13E, situated parallel to the rim of the Manson disturbed area, a 24 trace 
record shooting north yielded first, second, and third layer velocities of 
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1761, 3982, and 5360 m/sec. The second shot detonated between the cables 
produced records in both the north and south direction. For the north 
spread, layers of 2061 and 5417 m/sec were detected and for the south, 
layers of 1762 and 4540 were found. The west end of station 6-1 4E was 
adjacent to the south end of station 6-13E. Here, the west spread detected 
layers of 1595 and 4764 m/sec, and the east spread detected layers of 1679, 
3077, and 6275 m/sec. It may be seen that neither the second or third 
layer may be easily correlated with any other layer on any other record. 

The graph of the combined distributions from the first, second, and 
third layer reflects all the features mentioned previously. The few points 
added at higher velocities are from deeper refracting horizons in the 
carbonate sequence of bedrocks in the eastern undisturbed area. Both the 
disturbed and undisturbed velocity distributions are suspect of being 
quite scattered. 
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VELOCITY CONSIDERATIONS 

The velocity distributions are not a true representation of the 
velocities of the seismic refractors. The presence of dipping refractors, 
velocity anisotropy, velocity reversals, hidden layers, and inadequate 
velocity contrast will affect, in varying degrees, the measurement of 
seismic velocity. Each will be considered in this order. 

Dipping Refractors 

At nearly all locations the shooting pattern was such that travel 
times in opposite directions were recorded. In the split spread configura- 
tion, where the shot was detonated between the two geophone cables, the 
travel times represent samples of data collected from separate geologic 
settings, and in the reverse profile configuration, where shots were 
detonated near the far ends of the cables, the travel times represent 
samples of data collected from one geologic setting. The split spread 
configuration will then produce 2 seismic records based on 2 separate 
sections of material, and the 2 velocities, taken from the first line seg- 
ments in each record, may be compared. Similarly, the velocities, taken 
from the second line segments in each record, may be compared, and so on. 
The reverse profile configuration also allows these velocity comparisons to 
be made and further allows a check to be made on travel time reciprocity. 

These two tests must be applied to each pair of line segments from 
corresponding records before the dip calculations may be made. These tests 
have not, to the authors knowledge, been suitably framed in statistical 
analysis as they must, if the validity of the ultimate results are to be 
considered in the light of error of measurement. The statistical tests of 
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coefficients for several regression lines is a topic of analysis of 
covariance, and may be found in many books on linear regression, for 
example Wine (1964). 

First consideration must be given to the question whether the slopes, 
b 1 and b 2 , of each pair of corresponding line segments are statistically 
different. If there is no statistical difference between the slope 
estimates, they may be combined, and one slope value may be computed for 
that pair. If they are different, then a valid pair of slopes whose 
reciprocals are apparent dip velocities may be analyzed. 

Following Wine (1964,. p. 554), 2 true model equations, based on 
separate sets of time-distance data, may be written as 

t li = a l + e l d li + e li 
t 2i = a 2 + 0 2 d 2i + e 2 . 

The least-squares estimators of a-j , a 2 , are a-j , a 2 , and the estimators of 

3-j, $2 are b '|» bg. Now under the assumption that e-j. and are normally 

2 

and independently distributed with mean zero and constant variance and 
2 

a 2 at any point on each regression line, the estimated regression of time 


where i = 1 ,. .. ,n-j 

09 ) 

where i = 1 ,. . . ,n 2 


on distance is 
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t 2i = a 2 + b 2 d 2i 


( 20 ) 


because the expectation of the random component of each equation is zero. 

Let the terms e-j, e 2 be the measure of deviation at each d.. for each 
regression line in Equation 19 and are given as 


e li = t li ‘ a l " b l d li 
e 2i = t 2i " a 2 - b 2 d 2i 


(21) 


These are, in a sense, a measure of randomness at each d. for each regres- 
sion line. The residual sum of squares of Equations 21 are 

n l 

RSS, = 2 e?. 

1 i=l 11 


rss 2 = 


n 


2 2 
E e 


i = l 


2i 


( 22 ) 


The test of the equivalence of 2 slope parameters is stated as the 
hypothesis H q : ^ = B 2 = $o* * Jnc * er i^e nu ^ hypothesis that 3-|-$ 2 is 

normally distributed and the common population variance is unknown, the 
statistic t is given by 


x 


s [(1/SSd-j ) + (l/SSd 2 )] 1/2 


(23) 


where 
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ssd 1 = ^ (d^ - < dl >) : 


SSd 2 = s (d 2i - (d 2 » : 


(24) 


? 

and the pooled estimate of variance s in Equation 23 is 
2 RSS 1 + RSS 2 
S " n l + n 2 ' 4 


(25) 


Now instead of applying Equation 8 to segments of one regression analysis, 
let it be applied to these 2 regression lines from separate records. 


SSd 


1 = RSS 1 / [(n r 2) of (b)] 


SSd 2 = RSS 2 / t(n 2 -2) of (b) ] 


The statistic x is 


(26) 


T = 


RSS-j + RSS 2 
n l + n 2 " ^ 


'(n r 2) o-j^(b) (n 2 -2) o^b)" 1 


RSS 


RSS, 


W 


(27) 


The t is distributed with n-j + n 2 - 4 degrees of freedom. 

After t has been calculated it can be compared with the Student's t, 
here labeled t, with the same degree of freedom and an appropriate level of 
significance. If x is greater than Student's x the hypothesis is rejected, 
and it is concluded that b^ and b 2 are indeed different. 

If the hypothesis cannot be rejected, the pooled estimate of 6 is 


given as 
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av _ SSd-j - + SSd2 - b2 

SSd 1 + SSd 2 

and (a^) and <a 2 ) are 

(a-|) = (t-j) - (b) (d-j) 

( a 2) = (**2) " 


(28) 


(29) 


In this case the analysis may not proceed to the second test with this 
particular pair of line segments. 

As an example of an analysis of slope differences, consider the case 
where two records from a reverse profile configuration had been best fitted 
with 3 line segments. The 2 first line segments would be tested to see 
whether or not there was a significant difference between the slopes. Then 
the 2 second line segments would be tested to see whether or not there was 
a significant difference between these slopes. And finally, the 2 third 
line segments would be tested. 

A fundamental criteria in seismic refraction is that opposite ended 
shots and receivers may be interchanged without affecting the travel times, 
in other words, they are reciprocal. From this, each pair of line segments 
crossing the opposite-intercept must have the same time value. To reduce 
this to a statistically convenient question, the line segments from a 
hypothetical reversed profile are replotted on one set of axes. Then 
in an n-layer case, the first line segments should superimpose, and the 
second to the nth matched pair of line segments should each intersect at a 
distance S from the origin, which is equal to the shot point to shot 
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point distance. 

Now by shifting the time axis to a new time axis at distance S, the 
second to the nth matched line intersect on this axis. The regression 
equations then become 


<> 



t 

2i 


+ b^S + b i d-| 

+ b^S + b2d 


(30) 


By shifting axis, the data are simply coded, and no statistical changes have 
been made. The statistical test that the two regression lines intersect 
at this new axis is Hq.* a-| + b^S = a 2 + bgS. The statistic t is then 
stated as 

T.. ‘ - 1 - ' ** * S(bl ' b2) m (31) 

s [ ( Sdd-j/n-jSSd-j ) + (Sdd 2 /n 2 Sdd 2 ) ] x 

n l n ? 

2 ^2 

where Sdd. = Z d,. and Sdd 0 = Z d„. . As was done with Equation 8, 

1 i=l 11 1 i=l ^ 

Equation 11 may be rewritten in terms of 2 regression lines from separate 
records as 

Sdd^ = n^a-j 2 (a)/o-| 2 (b) 

(32) 

Sdd 2 = n2a2 2 (a)/0 2 2 (b) 

Then by substituting the Equations 32 into Equation 31, the statistical 

x is 

0 
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3-j - 



As before, the calculated t is compared with Student's t with n-j + n 2 - 4 
degrees of freedom at an appropriate level of significance. If the t is 
greater than Student's t, the hypothesis that the two regression lines 
intersect at distance D is rejected, and therefore, there is a significant 
difference in the opposite intercept times. If the hypothesis cannot be 
rejected, then there is no significant difference in the opposite intercept 
times. 

Because the data from split spread profiles were collected from 
different geologic settings, only the test of slope may be applied. And 
where reverse profile coverage is available, dip may be calculated on the 
top of the second layer only when (1) there is no significant difference in 
the slopes of the first layer, (2) there is a significant difference in the 
slopes of the second layer, and (3) there is not a significant difference 
in the opposite-intercept times of the second line segments. 

The line segment pairs for each seismic location were tested for slope 
and opposite-intercept criteria. Unfortunately when the least square fits 
were calculated only the combined residual sum of squares was retained for 
those line segments connected by a c-join. Therefore no estimate of the 
residual sum of squares for the left and right hand line segments was 
retained. Because of this, neither test of covariance could be applied to 
any combination of lines where one of the lines did not have a calculated 
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residual sum of squares. Also rejected from consideration were those line 
segments which had been fitted through only two points. 

The results of the slope test are shown on Figure 8. A significant 
difference in slope was inferred at the 90% confidence level, or simply 
that there was at least a 90% chance that there was a true difference in 
slope. Of the stations tested for slope difference, there were approxi- 
mately 30% (27) where there was a significant difference in the first 
velocity and approximately 70% (58) where there was not. The results have 
been drawn in the form of branches to show the different possible cases. 

Dip calculations on a 2-layer model.could only be made on those sites 
where there was not a sufficient difference in the slopes of the first 
layer, but where there was a sufficient difference in the second layer. 
Where no differences were found in either the first or second layer slopes, 
a horizontal interface model was assumed. This does not imply that the 
interface is truly horizontal, but that the magnitude of the subsurface 
undulations of the interface are of such a sufficiently low relief that the 
time deviations due to this surface have been combined with the error of 
measurement of time so that the model is simply a plane passing through the 
average depth to this surface. It is interesting to note that only 23 
stations passed the slope criteria for the 2- layer dipping model , and of 
these, only 1 case could be tested for dip of both the second and third 
layers. 

The results of the opposite-intercept test of the 23 sites that had 
passed both first and second layer slope tests is presented in Table 2. As 
in the slope test, a significant difference in opposite-intercept times 
was inferred as there was at least a 90% chance that there was a true 
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Significant difference in slopes 

First layer Second layer Third layer 

(b v b 2 ) Oy.b^) ib } \b z ") 



KD 


a Total of all stations. 

“Station with good borehole control. 

Figure 8. Slope test of all seismic sites by line segment 
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Table 2. Opposite-intercept test of 23 sites that had passed 2-layer slope 
test for a 2-layer dipping layer model 


Seismic 

site 

Shot 

no. 

Significant 

difference 

Relation to Manson 
disturbed area 

Outside 

Edge Inside 

2-2W 

1+2 

Yes 



3-3E 

3+4 

No 


X 

5-5E 

1+2 

No 


X 

5-1W 

2+4 

Yes 



5-6W 

6+5 

No. 

X 


6-1 1 E 

1E+1W 

D 

"h 



6-6E 

1+2 

D 

K 



6 -lE a 

6+5 

D 

k 



6-4W a 

1+2 

D 

u 



6-7W a 

3E+3W 

D 

u 



6-7W a 

3E+4 

D 

k 



7-6E 

4E+4W 

D 



7-2W 

4+3 

No 


X 

7-3W 

1+2 

Yes 



7-4W 

1+2 

No 


X 

8-2W 

5+6 

Yes 



8-5W 

2+1 

No 

X 


9-4E 

1+2 

No 


X 

9—1 W 

1+2 

Yes 



9-3W a 

6+5 

No 

X 


10-3'E 

3+4 

No 


X 

0-1 E 

7+8 

No 


X 

10-2W a 

9+10 

No 


X 

11-2W 

8+7 

No 

X 




12 

4 

5 3 


a Sites with well control. 
^Not calculated. 


difference in time. Only 3 seismic locations within the Manson disturbed 
area passed all criteria on which the 2-layer dipping refractor model is 
based. 
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Velocity Anisotropy 

Of the 46 reverse profile stations which were within the disturbed 
area, there were 27 stations at which there was no significant difference 
in the first layer; 12 stations where there was a significant difference 
in the first layer; and 7 stations where no determination could be made. 

At the 12 locations where there was a statistical difference in the 
first layer, there were 6 locations where there was a significant differ- 
ence in the second layer and 6 locations where there was not a significant 
difference. Significant velocity differences in the first layer are 
explained by the great variability of the unconsolidated drift. Lenticular 
sand bodies within the drift, variability of both grain size and lithology, 
and degree of water saturation may have an effect on glacial drift veloci- 
ties. 

At the other 27 stations where the slope test had been made, there were 
7 stations where tests of significance in the second layer could not be made 
and 12 stations where there was riot a significant difference in velocity. 

At the remaining 8 stations, the slope tests for dip calculation were 
successful. Finally, of these 8 stations, 5 stations did not pass the 
opposite-intercept test and 3 stations did pass. Only at these 3 sites may 
a simple 2-layer model with a dipping interface be applied. This is approxi 
mately 9% of the total number of stations within the Manson disturbed area. 
Those stations where significant velocity differences in the second layer 
were found, indicate that the first bedrock layer was anisotropic to seismic 
waves. 

Velocity Reversals and Hidden Layers ! 

Before the areal distribution of the first bedrock velocity may be 
considered, the presence of a layer between the glacial drift and the 
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bedrock must be considered. If the velocity of this material is lower than 
the velocity of the glacial drift, the layer will be totally undetected by 
the horizontal refraction method and will only be revealed by borehole 
seismic measurements. The effect of a velocity reversal is to reduce the 
velocity of the upper layer. In other words, an average velocity is 
obtained for the first 2 layers. Because the first layer for this area of 
study is glacial drift, the presence of an appreciable number of cases where 
a second layer with a velocity lower than that of the glacial drift would 
skew the distribution to the left. From the distributions of the first 
layer velocity in Figure 7, this is not the case, with the possible 
exception of one occurrence of a first layer between 1450 and 1500 m/sec. 

Green (1962) has investigated the case where a layer has a velocity 
between the velocity of an overlying layer and the velocity of an under- 
lying layer, but whose thickness is below a certain critical thickness. In 
such a case the layer will not be detected; therefore both the velocity and 
the thickness will not be measurable. Hence it is a hidden layer. There 
are also cases where a refractor is not detected due to the spacing of the 
geophones. In the former case, the refractor is a true hidden layer because 
it will never be recorded regardless of the geophone spacing. In the later 
case, the refractor is an apparent hidden layer. 

A true hidden layer between the first and second layer, whose velocity 
was intermediate between the two, was suspected at six locations (3-3E, 

4-1 E , 5- IE, 5-4E, 8-6E, and 11-1W) (Figure 1). Unfortunately well logs from 
boreholes near these sites did not reveal the presence of either a recog- 
nizable basal lithologic unit in the glacial drift or an unusual first 
bedrock unit. 
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At locations 3-3E, 4-1 E, and 8-6E a hidden layer velocity of 2300, 

2200, and 2250 m/sec, respectively, was assumed. These were chosen by 
averaging the second layer velocities from stations east and west of these 
sites. Calculations based on Green's method were performed on the veloci- 
ties and intercepts at these sites, and the maximum thickness of a true 
hidden layer for each record was determined. 

At sites 5- IE and 5-4E results from shots in one direction were 
inconsistent with the rest of the records and were discarded. At site 11-1W 
no hidden layer velocity was assigned because the velocity of the Pennsyl- 
vanian shales and sandstones in the area have not been measured accurately. 

Velocity Contrast 

Once a reasonably accurate representation of the first bedrock layer 
has been established, the areal distribution of the velocities will 
reflect, in some fashion, the lateral variations of the bedrock material. 
These results are presented on the left side of Figure 9. The graphs were 
computed by averaging the data by refraction line. There appears to be no 
trend in the averaged first layer (glacial drift) velocities. The averaged 
glacial drift velocity of 2143 m/sec on line 4 is unexplainably high. The 
second layer velocity for individual velocities below 3000 m/sec is pre- 
sented to illustrate the variations in the disturbed shale bedrock. Rather 
prominent regions of higher velocity are present on line 2 and 9. The 
utility of the algorithm developed previously is demonstrated on line 4 
where there is a lack of velocity contrast between the first and second layer 
velocities. In such cases the statistical determinations of the best fit 
line segments may be the only method by which consistent velocity data may 
be obtained. 



igure 9. Bedrock velocity graph and contour map of the Manson disturbed area 
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BEDROCK VELOCITY AND TOPOGRAPHY 

The mean velocity of the two first bedrock velocities, measured at 
each station, was chosen to approximate the true bedrock velocity. Here, 
the first line segment is chosen to represent the velocity of glacial drift, 
and the second line segment is chosen to represent the velocity of the first 
bedrock layer. This calculation is always higher than the true velocity of 
the bedrock, if all requirements for dip calculations have been met on the 
seismic plots. The equation of the true refractor velocity, as presented by 
Nettleton (1940, p. 269), may be rewritten in terms of mean velocity. Let 
V be the true refractor velocity, 9 be the apparent dip angle, and <V) be. 
the mean velocity of the refractor as measured from opposite directions at 
a seismic station, then V = [4R/(1+R) ]cos 0(V), where R is the ratio of 
the lower refractor layer velocity to the higher refractor layer velocity. 
The measurement (V) is less than 5% greater than V for R>0.68 and 9<16°. 

For all stations in the Manson area R>0.65. Therefore, the mean velocity 
of the first bedrock layer is less than 5% higher than the true first bed- 
rock velocity. 

The areal distribution of the first bedrock layer is shown on Figure 9. 
The area in the center of figure bounded by the 3500 m/sec contour roughly 
coincides with the area of crystalline rock, as determined by borehole 
methods. Another area of high bedrock velocity is located 10 km due west 
of the center of the Manson structure and has a diameter of 2 km. It is 
reasonable to assume that the crystalline rock also subcrops in this area. 
The limestones and dolomites which comprise the bedrock on the eastern 
flank of the undisturbed area are delineated quite well by the 4000 m/sec 
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contour, but velocity data alone may not be used to delineate the edge of 
the disturbed area elsewhere. Borehole information near location 11-1W 
reveals that the first bedrock layer, Pennsylvania shales and sandstones, 
is approximately 3 m thick in this area. The bedrock velocity at station 
11-1W reflects the velocity of the underlaying Mississippian limestone 
and sandstones, and the Pennsylvanian shales and sandstones have been 
averaged in with the glacial drift. 

The determination of the depth to the seismic refractors may be calcu- 
lated by a variety of methods. A method employing the concept of delay 
time, introduced by Gardner (1.939), has been successfully applied [for 
example the study by Pakiser and Black (1957)]. A similar method employing 
the concept known as time-depth was developed by Hagedoorn (1959) and. 
Hawkins (1961). Both of these methods make some assumption as to the ray- 
path. Because the bedrock is suspected of having velocity anisotropy and 
an irregular surface, depth determinations were made with the velocity-time 
intercept method, presented by Nettleton (1940), Dobrin (1960), and 
Steinhart and Meyer (1961). Because very few seismic records could be 
subjected to dip calculations, depth calculations were based on a simple 
horizontal model. 

In view of the obvious simplicity of this model, error bounds on the 
depth determinations were made purposefully high. Let the first layer 
correspond to the first line segment, j=l. The error bounds of depth to the 
top of say the jth^ layer (j>l) were calculated with error bounds on the 
parameters of the jth^ line segment. Confidence limits were chosen so that 
there was a .99 probability that the true time intercept and wave slowness 
were within the calculated bounds. The time intercept and wave slowness 
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for the first to the jth^ line segment were first used to calculate the 
depth to the top of the jth^ layer. The lower depth bound was found by 
substituting the value of the upper confidence limit of the time intercept 
and the value of the lower confidence limit of the wave slowness of the jth 
regression line and then by recalculating the depth. The upper depth bound 
was found by substituting the value of the lower confidence limit of the 
time intercept and the value of the upper confidence limit of the wave 
slowness and then by recalculating the depth. The upper and lower equations 
for the jth_ line segment with n data points are 


tl u = a. + TOj(a) * tbj ♦' Wjfbjjd, 

‘i 1 = a j ' TO j (a) + [b j + TO j< b > ld i 


(34) 


where the superscripts u and 1 indicate the times for the upper and lower 
line segments and where t is Student's x Q with (n-2) degrees of freedom. 

Figure 10 is based on an interpretation of a combination of the areal 
distribution of the depth determinations and bedrock elevations from bore- 
holes. The irregular surface within the Manson disturbed area is sharply 
contrasted by the gentle slopes outside the area. The average elevation, 
within the boundaries of the disturbed area is about 300 m above sea level, 
whereas the 320 m contour line roughly traces the average elevation of the 
undisturbed bedrock surface outside the boundaries. Within the disturbed 
area a number of bedrock highs and lows are prominent. The central area 
of crystalline bedrock is seen to be a topographic high. Hills of disturbed 
shale in the northeast and southwest portion of the disturbed area are 
present. Six closed depressions surround these bedrock highs. The 
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drainage ways of the disturbed area are not well developed. Several 
undisturbed bedrock valleys terminate in the depressed area. Two bedrock 
valleys in the southern portion trend northwest-southeast and their gradient 
is toward the southeast, away from the disturbed area. The eastern bedrock 
valley is a headward extension of the bedrock valley shown by Hale (1955, 
p. 151) to extend southeast across the southwest corner of Webster county. 
The major drainage network in north-central Iowa is known to drain to the 
southeast. 
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DISCUSSION 

A number of features have come to light from the seismic investiga- 
tion. The Manson disturbed area has been shown to have the appearance of 
an impact crater. A central uplift of extremely deformed crystalline rock, 
which was lifted some 700 m from the basement, now extends through a 
shallow circular basin of highly deformed Cretaceous shales. At the time 
of burial by Pleistocene drift, the area surrounding the uplift had 
not developed an integrated drainage network, but occurred as a number of 
rises and shallow closed depressions. The lack of an integrated drainage 
network within the disturbed area may indicate that the structure is 
relatively young. 

The relative erosional resistance of the Mississippian limestones and 
dolomites beyond the northeast rim of the area contrasts with the lack of 
resistance to erosion of the Pennsylvanian shales and sandstones beyond the 
southwest rim. This is expressed by the very steep slope along the rim in 
the northeast and the more gentle slope along the south and southwest rims. 
Two rather extensive headward erosion channels had developed in the Pennsyl- 
vanian shales and sandstones prior to burial. The western bedrock valley 
follows the edge of the disturbed area for some distance. The zone between 
the disturbed and undisturbed bedrock must have been less resistant to 
erosion than either the disturbed or undisturbed bedrock. The shape of the 
western bedrock valley suggests that the disturbed bedrock was comparatively 
more erosionally resistant than the undisturbed Pennsylvanian shales. 

The second crystalline subcrop, as suggested by the bedrock velocity 
map, is also seen to be a western extension of the center crystalline 
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bedrock high. Indeed, the gravity high (Figure 3) in the area coincides 
with this suspected subcrop of crystalline rock. The elongate 4605 gatima 
magnetic anomaly (Figure 2) is situated parallel to the crystalline bedrock 
high and shifted 1 km to the east. 

Only 20% of the seismic records taken at sites well within the Manson 
disturbed area could be represented by a 3-layer seismic model. At many 
of these sites the third layer detected by a shot in one direction was not 
detected in the shot record in the other direction. A third layer with 
velocity about 3400 m/sec may be present below sites 3-2E and 3-1 E, but the 
depths to this layer are quite irregular and dip calculations could not be 
made. For the most part, there is no mappable seismic refractor, below 
that of the drift-shale interface, to a depth of 250 m. 

The rim structure appears to have considerable complexity. The line 
segments fitted to the records at stations 6-13E to 6-16E showed that there 
is a large variation of velocity in the undisturbed bedrock east of the 
Manson disturbed area (seen earlier). No consistent structural model, no 
matter how general, could be made of the geometric arrangement of the bed- 
rock layers. Simple seismic models may not be applied in these areas, and 
it points to considerable geologic complexity, probably due to complex 
folding of the bedrock units adjacent to. the Manson disturbed area. 

Another significant feature is the seismic velocity of the crystalline 
rock in the center and western portion of the structure. Press (1966) has 
reported that granite velocities may range from 4.8to 5.9 m/sec. Velocity 
histograms of granite and gneiss were presented by Hasselstrom (1969) for 
shallow seismic refraction measurements. In each case at least 75% of the 
sample had a velocity greater than 5000 m/sec. At Manson the velocity of 
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the mixture of granite, gneiss, and diabase is between 3000 and 5000 m/sec 
and averages near 4000 m/sec. The unusually low velocity of this crystal- 
line mass may indicate that this rock body has an abnormally high porosity. 
This may be a result of the large scale brecciation of the crystalline mass 
during the implacement event. 

The unusual deformation characteristics of the minerals of the central 
crystalline mass have shown that this structure could have been created by 
one of only two 'catastrophic' events. Several similarities of the Manson 
structure to the Canadian craters are noted. This crystalline rock has a 
close similarity to rocks taken from central uplifts in larger Canadian 
craters, and the age of this crystalline mass is near the approximate age 
of the basement rock as observed in some Canadian craters (Michael Dence, 
personal communication, 1970). 

Regardless of which mechanism had created the structure, the basement 
rock would be involved in the final disposition of the materials, because 
the thickness of the sedimentary column is considerably less than the 
diameter of the structure (roughly 1:400). Therefore gravity and magnetic 
anomalies from the upthrown basement rock may not be used to distinguish 
between the two mechanisms. 

The questions raised by the occurrence of this complex structure must 
be answered by sifting and evaluating a wealth of data, some of which is 
contradictory or coincidental. The nature of the origin of the Manson 
disturbed area must await the collection and interpretation of additional 
data. 
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APPENDIX A. LIST OF THE COMPUTER PROGRAM 

1 . Purpose ; 

The program fits joined regression lines to an ordered series of 
measurements on the basis of residual sum of squares. For each line seg- 
ment, the parameters and their standard deviations are calculated. The 
measurements are taken to be seismic refraction time-distance pairs, and 
from the results of the line fitting analysis, depths and depth error 
bounds are determined. 

2. Method ; 

The set of up to 24 time-distance measurements are sequentially 
partitioned into subsets of data, each of which is searched for the best 2 
segment fit. Each subset has an upper bound, which is chosen prior to the 
analysis, and a lower bound, which is determined by the position of the 
join point from a preceding 2 segment fit. The final solution is derived 
from a combination of the 2 segment fits. 

For a subset of data, pairs of regression lines are fitted to all 
combinations of the data, and the pair of regression lines which have the 
least combined residual sum of squares are taken to be the best fit. 

Least squares regression is used to best fit the line segments. The 
intercept and wave slowness (parameters a and b) of each regression line 
are calculated, and under the assumption that the variable, time, is inde- 
pendently and normally distributed with constant variance, the standard 
deviations of intercept and wave slowness are found. The depth calculations 
are based on the standard time intercept-wave slowness equations. Error 
bounds in depth have been treated previously. 
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3. Accuracy : 

The accuracy of the program was measured from test data presented in 
Table 3. These data were analyzed by (1) the OMNITAB program available at 
the Iowa State Computation Center, (2) a least squares program written by 
Bevington (1969), and (3) a subroutine from this program, SJOIN. All 
programs were run on the Computation Center's IBM 360/65. An evaluation of 
widely used regression programs by Wampler (1970) found that the single 
precision OMNITAB program is correct to 4 and 6 significant figures for the 
first two coefficients of a fifth degree polynomial. From the test data, 
the calculations of a and b by subroutine SJOIN agree with those from the 
OMNITAB program to 5 and 6 figures (Table 4). The standard deviations of 
a and b agree to 7 and 6 figures. The calculations of residual sum of 
Squares from the subroutine SJOIN and program OMNITAB are included. The 
results of the Bevington program are presented as another comparison. 

Because the calculation of residual sum of squares is fundamental to 
the validity of the program, the time values, t (Table 3), were adjusted 
for a range of values. The results are presented in Table 5. The approxi- 
mate range of values used in the analysis of the data from this study were 
0<d<800 and 0.005<t<0.450, so the residual sum of squares calculations are 
valid to approximately 5 significant figures. Because any 2 calculations 
of the residual sum of squares were equal to at most 4 digits, the program 
is considered to be valid for the statistical decisions undertaken. 
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Table 3. Test data 

of distance, d, and time, t 

d 

t 

15. 

12. 

90. 

28. 

195. 

45. 

300. 

63. 

405. 

82. 

Table 4. Comparison of three computer programs 


Coefficient 


Program 


SJOIN 9 

0MNITAB b 

Bevington 9 

a c 

10.53826 

10.538361 

10.538420 

o( a) 

0.91706004 

0.91706008 

0.916853 

b 

0.17642574 

0.17642599 

0.176426 

a(b) x 10“ 2 

0.37415442 

0.37415458 

0.3741 

RSS d 

4.1312900 

4.1345215 

-- 


a Double precision. 

^Single precision. 

c The form of the equation is t = a + bd. 
^Residual sum of squares. 
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Table 5. Significant figures of the calculation of residual sum of squares 
for various values of time 


Data 

Approximate 
veloci ty a 

Significant 
f i gures b 

<W C 

5.833 

>13 

(d- ,0.1t- ) 

58.33 

5 

(d.,0.01t.) 

583.3 

6 

(d. ,0.001 t.) 

5833. 

4 


a Velocity is the reciprocal of b. 

b Based on the manual calculation RSS = 4.131290027447. 

Slhere i = 1 ,. . . ,5. 


4. Usage ; 

a. Input ; 

A standardized punched card record format is described in the 
conment cards of subroutine RREC. 

b. Output ; 

Although the length of the output is a function of the number of 
fitted line segments, the number of data, and the number of acceptable 
regression candidates, a typical 2 segment analysis of 24 data points will 
yield about 1 page of printed output. 

c. Operating Procedures : 

The first card of the data deck must specify the mode of operation 
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of the program (see the comments in the main program). The next 2 cards 
must be the first 24 Student's values with a format of D6.3. The input 
records then follow. The end of data is specified by placing a blank card 
and a card with 99 punched in cols. 47-48 after the final record. 

d. Constants : 

A total of 4560 bytes has been allocated to unlabeled COMMON and 
116 bytes to COMMON area /LABEL/. 

5 . Subroutines Requi red : 

The subroutines required are RREC, SETS, SJOIN, CJOIN, OUTPUT, and 
ABYSS. The main program calls subroutines RREC, SETS, and ABYSS. Sub- 
routines SETS calls subroutines SJOIN, CJOIN, and OUTPUT. 

A list of the program follows. 


71 


r-4<Mco^inOf'-coa'O*-<<Mcn^tn*or^coo'Or-«cyjfOtf-in%of^cD0'O«-«cMco<4-inNOi^ 

oooooooOOHHHHHHHHHH(\j(\jN(\JNMM(MwrorO(nnmi^^foro 

ooooooooooooooooooooooooooooooooooooo 


• 

• 

_J 

UJ 

o 

to 

o 

z 

z 

O 
•— » 

a: 

1- 

UJ 

a 

>- 

O 

<r 

o 

1 

z 

CM 




P -4 






*— • 






f- 










z 






z 0 






z 










< 






z z 






< • 










♦* 






<r < 



• 



















a: • 



O 



00 










%o 






O CC cc 



UJ 



z 








00 








000 



H- 



0 








O 


>0 






a: K k ♦ 



< 

• 


0 








£* 


w 






a 0 0 > 



_J 

0 










O 


<c 






UJ UJ J 



3 

UJ 


h- 








O 


K 






-J > > z 


• 

O 

1 - 










UJ 


-J 






< 0 



-J 

<r 


00 








cc 


UJ 






a: H K 



< 

-j 


• 










0 






U)ZZ J 



O 

3 


H- 








0 


*«• 






>HWUJ 




O 


Z 








UJ 








UJ O O O 



K 

-J 


UJ 








0 


CM 






00 Q- Cl O 

• 


O 

< 


0 








3 

to 

CM 






Z 

0 


Z 

O 


3 








O 

00 



CO 




X h h 

UJ 





K 








UJ 

>• 

in 






K to 00 oc 

X 


UJ 

UJ 


CO 








cc 

QD 



K 




H UJ tu UJ 0 

• 

oc 

CC 











< 

< 


Z 




3 h h > 

z 0 < 

< 


UJ 








0 

•» 

z 


HH 




<1 

3 UJ 




X 








K 

to 



a: 




cl a: ol j a 1 

to 

CO 


H- 









K 

< 


CL 




UJ UJ UJ | 


O Z 

z 






A 




</> 

ID 

cd 


z 




to a a cm k- 

z 

0 

0 


UJ 




3 




to 

00 

•» 


•» 




3 

O 3 

*— 1 

•-H 


CC 




1-^ 




>- 

•» 

*»• 


v. 




00 00 00 z a 

K 

1 - 


< 

< 4 - 



#• 




CD 

0 

ro 


CM 




UJ < <c < 



< 

<t 



CM 



O 




< 

UJ 

CM 


*V 



to 

X 

00 to Z 

z 


00 

•» 



H 





cc 

•* 


X 



K 

*- 0 0 0 

O O «-H 

1-4 


a 

H 



•» 




O 

oc 

00 


0 



CL 

UJ UJ UJ 

CC 

OC 

z 

z 


oc 




O 




Z 


w 


z 



z 

00 Kl fM M 

<r 

< 

oc 

OC 

z 

«r 

II 



•» 




< 


0 


3 



*■ 

UJ > >- > 

0 CJ 

ID 

UJ 

sc 

u 




K 





•• 

*- 


CL 



— «< 

O -J -J -J 



K 

K 

•» 


HH 



•> 

O' 



00 

0 

CO CM Z 



z 

*-* < <t < 

> > UJ 

UJ 

-J 

0 

•- 



O 

O' 



h- 

ID 



► 



CM Z 

> z z z 

H- 

K 

O 

0 

*c 

CL 




Z 




UJ 

CC 


to V. 



< * 

0 < < < 


*— « 



•* 

►-« 

1 



UJ 0 


0 

00 

*-» 

CO 

•* 

v-i 



M (j) ^ 

CL 

O O 

X 

X 

* 

X 

w 




H 


H 


3 

(M 




a 

•« •* 

a < < <t 

O O H 

K 

*: 

H- 

CO 



•» 




> 

O 

*- 

in 0 


O' 

- O O 

H K K 

-J 

-J 

CL 

CL 



w 



-J O 


«H 

_J 

UJ 

CM 

•» 

< 

in * 4 - 

H> HH 

a < < < 

UJ ID 

UJ 

UJ 


O 

CO 



-i CD 



a 

a: 

H sO UJ 

w 

w- 

OlOM 

0COOQ 

> > O 

O 

00 

Z 

A • 

CM 


«• 



11 

a 


w- 

w 

a: 

a 0 

H •» _J 

< 





•* A 

< 

0 0 



-J 




< 

00 

< < z 

♦-H 

•» 

w < UJ 

O 





0 1-* 


- 0 

H 



O' 


K 


tu 


H 




DhCD 

0 

1 

2 
0 

f-C O 


c ^ 

O 

O CM 




O' 


00 

0 

z 

UJ UJ >J* 

** 

O' 

mUJ < 

f— 





uj ro 

z 

< r-4 

II 

P-4 

u 

• 


UJ 

H- 


CD 

M * 

* 

-*• 

CD — 1 

00 II II II 

II 

II 

II 

II 

oc ^ 

0 

UJ w- 



UJ 0 


H 


K 



oc 

CC 


a: \ 

CC 





z 

0 

OC 

l-H 

II 

a' 

UJ 


•-H 

CL 

3 

00 

00 

UJ 

UJ K 

UJ z z 

*-* r-4 

CM 


CO 


w h- 

UJ 

Z K 



0 L 

• 

0 


3 

O 

* 

* 

0 

CD 


CD O O 

U. 





< 

00 

w < 

O 



-1 


0 


CC 

— J 

-J 

UJ 

UJ 

-J 

UJ Z Z 

• 

• 


• 


0 z 


O Z 

f-H 


-J 

w 

II 

CM 

H 

CO 

< < 

K 

1 — 

<: 

K Z Z 

UJ -J 

-J 


-J 


< oc 

UJ 

< oc 



-J 




UJ 

3 

UJ 

UJ Z 

Z UJ 

ZOO 

X 0 

O 


0 


UJ 0 

X 

UJ 0 

0 •-> < 

UL < O 

00 

00 

cc 

cc 

H-t 

h-i a: 

MUO 

1 - 0 

O 


0 


a: u. 

1 “ 

cc u. 

0 ** 0 

P-4 

sc 0 














CO 


O' 


O 





u u u u o 


oououooooo ooo 


o 


72 


oo(7'0*-HC\irn^iri'0[^co(^o«^{\jm>a-mNor-a30>o«-H<\Jco>fru^NOr-‘COO'0*-HC\icO'4- 

rnm>t't'4 , '4’'t^'4*>t^>tiAirMniAir\iniAifMninsO'0>0'0>0'0'0'00'Oi w ^NNN 

ooooooooooooooooooooooooooooooooooooo 









0 

X 

h- 

0 

• V 
















H • 

CM 















2 00 

1-4 















04 h- 

•k 















CL 

X 















LL 

r- 















O 2 
















* 

in 







l 








O 0 

• 















I- X 

CM 















00 >o 

H 
















o 














Mi 

«k — 

0 

<— 













2 

X 00 

X 

00 













►-4 

4* OO 

M CM 

• 













2 

0 

M «k 

o 













< 

• Q 

w o 

CM 


«M 











0 

1- 04 

00 • 

o 


JD 











00 

a oo 

H in 

► 


0 ' 











I- 

UJ UJ 

Q- ll 

m 


H 











o. 

o oc 

2 














2 

CC o 

•* X 

II 


II 

-j 










0 

UJ 0 

—» in 




> 










<C 

H X 

4> •• 

to 


-o 

-j 











2 >t 

0 A 

UJ 


0 


O 









-J 

04 0 

r-l in 

CC 


<M 

UJ 










UJ 

m •» 

• 

< 


04 

X 










o 

0 >• 

II CM 

3 


0 

H 

0 


Ml 







0 

x w- 

H 

a 


-o 


K 


< 







-J 

►-! 

-o O 

to 



K 



►— 







UJ 

0 OJ 

0 



< 

3 

O — 


UJ 







O 

- O 

— X 

LL 


H- 

O 

O • 


hJ mi 







1—4 

CC -J 

i-4i CM 

O 


-J 

X 

II 


0 O 



- 




► 

UJ UJ 

•k W 



UJ 

h 

00 « 


< CO 







o 

> > 

-? CO 

2 


O 

04 

Mi W 


H- •* 

CO 






•» 

< 0 

w 0 

3 


w 


O CM 


-J H* 







H 

-J * 

< 4* 

00 


0 


• CO 


UJ 2 

o 






0 

X 

1- • 



-J — 

O 

Of r-l 


O H4 

H- 






3 

- 

-J o 

-J 


UJ 0-4 4* 

2 

UJ 0 


0 X 







1-4 

^ 0 

UJ *H 

<t 


o — in 

< 

• 0 


00 a 

o 


rH 




•» 

m 

QU 2 

3 


04 00 • 


2 


0 2 

o 






►** 

00 

— ' 0 04 Q 


CO 0 H- CM 

X 

X 


00 «■* 



II 




00 

<M 

- X 2 

0-4 


cn o a •-! 

H* 

• — N* 


fr— 

Ml 






UJ 

• LL 

I-H CO < 00 

< 

•-4 2 a 

04 

• *m 


CL UJ 

o 


H* 




h* 

* O 

0 

UJ 

H 

O 40 ^ 

X 

OC 'V 


2 h- 

0 


00 




•-4 — • 

0 

~ (\J — 

CC 

< 

H 0 CM * 


0 — • 


0 0-4 

a 


UJ 


Mi 


► O' 

W Q 

JHh w 


Q 

O w JHX 

X 

• 4* 1 


-J QC 

UJ 


l~ 


in 


O NNH 

UJ CO •* CO -J 


O in 0 uj M) cm 

H* 

m ro 0 


LU 3 

• 


04 


•* 


Z - fO </> — 

O 0 X ► 

< 

oc 

O 0 XO 00 

O 

r-l 0 W 


O 

«M 




»- 


UJ h- 

H • ^ 

h h CO h h 

UJ 

X CM H x CM 

GO 

• H CO 


HW0 

K 


— * 


2 


0H 2 

0 •» • 

0 2 -20 

> 

«— O 0 0 O w 


a 2 v 


00 CM 

00 

CM O 

0-4 

• 

W 1—4 

• X H H • 4-4 

»- 

< 

O 2 4* r* 2 * 

CO 

UJ 04 0 


+ • 

UJ t-J 


• 

Q al h 

(i O h- • 

OC QC 

O 

-J 

• 3 < 3 X 

X 

• OC r*4 


00 O 

h H 

+ 

a 

0 

Cl 

• 

00 CL 

m *k W 

II CL • CL 

0 


GP CL in II CL <* 

H- 

j a • 


X 00 UJ 

1-4 LD 


UJ o z 

w 

1- 2 

• <M 

Z40Z 

w 

a 

UJ 2 ~ 2 — 

a 

UJ 2 w 


X >- • 

^ 3 II 

< 

• 


— * 


UJ ^ 

UJ CO 

0-4 W w 


3 

• w 0-4 

UJ 

O w 

0 

X GO X 

3 2 

X X 

II 


•“ 00 

hQ-H 

h 

H* 


-J H- H- 

0 

0-4 H 


< X 

04 04 (— 


x 


UJ < 

UJ < o ► O UJ < LU < 

X 

X UJ < O UJ < 


w UJ < 

II 

II X 

w H 00 

II 

w 

2 h- 

2 

-J K 

2 -J • 

fOhZh 

2 

o 

w K 3 l 0 H- 2 

a 

w H 2 



2 UJ 



04 

»-0 

-1 ^ 

QC OO 

M 0£ H 

CC 

. 2 

0-4 OC 04 

2 

0M CC 

* X -i 

HOh 

<r 

LL 

X 

a: 

o < oc 

o - - 

O QC o oc 

o 

3 

LL QC O O QC O 

0-4 

U.0CO^^<ll 

w o w 

x 

1—4 

< 

3 u. o ^ 

LL 

OXU. 2 

LL 

a. 

wJILQJU. 

LL 

04 2 U. 

^ 0 •-< 








X x 









O 

ro 




in 


O' 

O r-l 

CM 


O O rH 


CO >t 


in 

CM 







CM 

CO CO 

CO 


in 0 


CO CO 


CO 


u o u 


u o u 


73 


irvNor^coO'0*-4Csjm^in'Or^ooo'Of-ifMm>tin>ci^coo'0»-*csJro^-inor^oDa'0»-4 

r*-r^i^r^r^coooooooooooaooococoo'0'0'0'0'0'0'0'0'Oo©ooooooooMiMi 

OOOOOOOOOOOOOOOOOOOOOOOOOm^mImImImImImImImImImImI 


V 


X 







00 


>4* 







H— 


M . 







CL * 


o 







X 


K 







Z 4- 


Mi 







o •* 


£ 

•* 

•* 





m m 


M- 1 

(NJ 

CM 





V. 00 


-J 

• 

• 





CM 



r- 

r- 





Mi U. 


cc 

UL 

UL 





- O 


UJ 

m 

m 








X 

X 





o 


o 

s0 






00 h- 


-i 

•» 

•» 





0C CO 



(NJ 

CNI 


CO 



UJ » 


•» 

• 

• 


O 



> * 

M* 

X 

I'- 

h- 


m* 



< X 

*0 

04 

ll. 

UL 


H 



_j m 

M 

m «“* 

•* 

m 


00 




H 

* • 

X 

X 


Mi 



UL • 


h- H- 

in 

in 


0C 



O UJ 

II 

1— * Mi 

•* 

•» 


LU 



a 


X X 

CsJ 

(NJ 


J— 



x o 

-3 

1— I Mi 

• 

• 


O 



UJ _J 

•» 

-J -J 

r- 

— 


< 



CD 00 



— u. 

O UL 


OC 



z . 

M 

oC X 

o ► 

* *> 


< 



— 3 m 

M 

LU UJ 

► X 

mI X 


X 



• Z X 

"3 

Q- D. 

Ml ** 

>4* 


a 



a: ^ 

w 

CL CL 

•» 

11 •• 





LU » •» 

< 

3 3 

II <\j 

O (M 


UJ 


M 

> • 

h- 

• « 

* 

Q Mi • 


> 



< st h- 

-J 

•» 

3 r*. 

co * r- 


< 


1 

-J < z 

UJ 

X m 

m LL 

o ~ il 


2 



in i-4 

a 

r\j x 

m» o» 

oo -o •* 


i 


< 

>- ► 

w 

m (NJ 

Mi X 

(Nl •» X 


CC 



K • LL — 

•* 

• * 

- o 

• M \0 


Mi 

Mi 

•* 

Mi _J O • 

M* 

to • 

3 * 

CO W * 


< 



O UJ to 

Mi — * 

00 K 

w cm 

< CM 



1 


O Q O O 2E 

w in 

UJ Mi 

< • 

* h- • 


UJ 


< 

JWQi f- 

</) • 

Z X 

h 

UJ h- 

o 

> 

< 

K 

UJ ►OlOQ 

K 04 

X Mi 

UJ UL 

— * rvl UL 

— r- 

< 

Mi 

-J 

> O U » »“» 

CL mI 

O -J 

ixj m 

-o ~ «* 


X 

w 

uu 

Mi LU * 00 

z o 

1— 1 

^ X 

•> X 

• o 


00 

o 

2 0C X LU 

o 

X a: 

in 

*-» in 

II K 

a 

H* 


O — • CO DC 

M» M -K» 

I- UJ 


W ** M 

» » 

H 

CL 

II 

J h ► ► • 

00 X O' 

• 2 

mI (Nl • O' 

< ^ (NJ 

in w o 


z 


co x - •* _i 

nMrn 

- o 

>t • oo co 

H* <4* • — 

sJ-MO 

a£ 

UJ 


X •* O H x LU 

X -J 

U. CC * LU 

uj * r- • 

•» CO 

UJ 

a ti 

Mi 


h h H fl. ^ O h 

(NJ 1- O' 

• 

H 

Mi J— LL 

UJ H 

h 'O N h 

U.hh 

Ml — 

> 

Ml (NJ Mi 

o + 

w Z *- LU •* 

M« Z 

H Z ► 

M 


H* Z - 

h- Z 

* * 2 

►UJ Z 

* OJ 

< 



•» M 

1 

2 

Hi » O • 

m Mi 

H M » 

X 

1 

Ml Mi X 

UU Mi 

»-l H II Mi 

X UJ Mi 

• * 

-J 

+ 

4 

Ml Mi 

Ml «J O 


CC O DC > H cC 

x a: o 


Q£ I-I X a. 

QC 

Ml UL QC O O 





CO uu 

• 

a. • uu h- 

CL 

Nt CL * 

•» 

-J 

II CLH 

CL 

ii it — a. 

-4 CL 

• LU 

H 

-1 

_) 

M 1 

II w Q || 


2 w H- ^ 

It Z 

^z^ 

» 

UJ 

Z M 

2T 2 

3 Z 

mZZ 

w • 

00 

UJ 

uu 

M Mi 


w z o 

W 


X Q Mi w 

Mi w 

O w * w 

Mi W 

X 

a: 

QQwcD 

-9 3 — 

K 

h* h O w 

K H* 

h- 

Mi 

K 

m 

Mi 

h • % 

K X 

Mi 

Mi 

Mi 

Mi 

~ II Ml 

< 

UJ < - _J 

UU 

< UJ < 

CL 


o UJ < 

•* UL) 

CO (M w LU 

< ► UJ 

< X 

UL 



O — 

O < w 

51 

1- X ► UJ 

>0 K 

I h Z UJ 

II 

^ H X 

X K 

^ >4* <t H 

ZXH 

X w 


II 

II 

<t (/) J 1 K J </) 

QC 

H O' X > 

M 

0C Mi QC 

a 


m ac o mh 

h Mi 

CC O w 

QC 

H- 



h- 

J UJ h 

O 0C O 4“ - 

o oc 

O CC O 

* 

LL 

O CC O Ml QC 

O O UJ oC 

O Ml QC 

O UL 

UJ 

4fflOa.OUJQb 

UL 

2 UL 

a 2 

UL 2 UL 


Mi 

O 2 LL 

2 

O O M 2 

u. 2 

U. «-h 

00 

Mi 

Mi 

Q Z O O *-* Z 


X X 



X 



X 


X 







O 

h- 

>o 

00 O' 



O .H 


04 CO 

^4* 

in 






(O 

m 


CO CO 





>4* * 

st 

-4- 





>4“ 


o o u 



74 


(NjfO>tu^vor^coa'0^(Mm>r^'Of^ 



O 

ULi 

X 

0 
z 
z> 

01 

z 

LU 

Ui 

CO 



CO 

UJ 



> 


o 

< 


H 

X 


O 

— 00 


o 

o a 



O CC 


*•» 

rH < 





rH 

V- 


• 

z >■ 

m 

in a> 


o 

r- uj 

a: *-* 

l 

• 

a o 

o 

o < 

z o 

in 

H X 

— -j 

o 


LU 

h- 

O • 

UJ > 

O O I-H o o o 

o o 

H- • 

o o o o o o 

z 

►H •* 

• • • » • n 

— < 

a: x 

O O CO o o o 

rH * 

x o 


t — * 

rH 

M II II II II II 

O rH 

a ** 


LU • 

rH • 


• a 

• 1 

H H H H s—4 H 

h- rH UJ 

Of m 


l/> • 

UJ ^ 

H M m ^ in nO fO 

LU || * 

rH • 


K ^ 

-J K 

<<<<<< o 

HH- 

O X < 


w CO w 

h- ~ x a 

— 1 » 1 1 1 — l 

UJ 

a: O 

OI LU ill UJ UJ UJ o 

li K tl 

OllOK 

O Q O O O O O 


o H LL 00 


o 

z 

LU 


o in o 

r- r- O' o 


u 


75 


H 

* * 

* 

* * 

* * * * 

* * * * * 

* * * 

* 

* * * 

* 

* 

* 

* * 

* 

* 

* * * 4fr 

* 

* 

* 

O 

* 


















O 

* 

1 


00 















u 

* 

O 


-J 















UJ 

♦ 

< 


< 









• 






a: 

* 

LL 


z 















a: 

* 

UJ 


CD 
















* 

a: 


»—« 
















* 

a. 


oo 















\ 

* 



















* 

Q 


H 
















# 

od 


•— i 
















* 

< 


















* 

o 


CK 






• 










* 



o 






oc 










« 

Q 


< 






o 










* 

2 


UJ 






H 










* 

< 

t 

a: 






O 










* 

x 

_} 







UJ 










* 

X 

— i 

oo 






> 










« 

O 

LU 

»-hI 
















# 

o 

2 







H 










* 


oo 

O' 






Z 










* 

LU 

O' 






>-« 










* 

X 

< 


a: 





O 










* 

H 


Z 

O 





Cl 










* 


00 

LU 

H 















* 


H 

X • 

O 





H 










* 

• 

2 

X H 

UJ 





oo 










* 

lo 


LU 

> 





UJ 










* 

a 

O 

* to 




UJ 


H 





O 



ad 



a: 

a 

UJ LU 

H 



o 







UJ «J 



O 


* 

o 


O < Q 

Z 



o 


ad 


• 



O UJ 



H 


* 

o 

H 

OHO 

•— * 



o 


o 


00 



O *-< 



O 


* 

LU 

00 

O < o o 

O ad 



UJ 


LL 


UJ 



O U. 



UJ 


*# 

0d 

UJ 

O 

CL O 



o z 




3 





> 

—* 

* 


H 

H > O 

OC H 



O O 


00 


-J ■ 



>■ H 




3 

* 

O 


UJ UJ cX 2 

O H O 



o +* 


LU 


< 



z ad uj 



LU 

►— » 

* 

•— 4 

O 

00 X H < 

H 00 LU 



H 


3 


> 



O H oo 



X 

•» 

* 

X 

H 

H 1U 

UUi> 



CD O 


-J 





t-H LU 



t-4 


* 

00 


< X H 

OC UJ H 


1 

Z UJ 


< 


O 



H X < 



H 


* 

HH 

o 

H LL O 

O > H 



►-« -J 


> 


9-1 



< O H 

1 



-J 

* 

LU 

H 

< O LU LL H LL Z 



X UJ 







U LU < 



UL 

*• 

* 

00 


O CD O O LU O h 



O 00 


LL 


LL 



o o 



O 

o 

* 


a 

O 

LU O O 



H 


o 


o 



LL 




•» 

* 

o 

3 

UL Z O X 

> Z X CL 



< II 







wQlL 



X 

H 

* 

LU 


OUJ<h 

< H 

a: 


CO 


co 





H < O 



H 

•* 

* 

O 

a: 

tu o 

LU H CD H 

< 


-J 


Od 


< 



Z UJ 



O 

Q 

# 

3 

o 

OUJoiZ 

X 00 Z 00 

X 


II UJ 


•-* 


H 



UJ OC o 



z 

Z 

* 

o 

LL 

zxauj 

M M UJ UJ 

UJ 


00 


< 


O 



a a z 



UJ 

LU 

* 

LU 


lUhWJhQ Jh 

0d 


X 1*1 


a 


H 


a i-4 oo uj 



-J 

HH 

* 

OC 

a 











CC 





•» 

* 


< 









LH 


< 





-J 

* 

00 

LU 



H 






CVJ 


o 





-J 

# 

o 

OC 



< 












•» 

* 

< 




X 






Z 


z 





-J 

# 

LU 

00 



a: 






X 


o 






* 

oc 

HH 



o 


a 




3 


►1 






* 





LL 


oc 




-J 


H 





u 

* 

o O 





<. 




O 


< 



00 


LU 

* 

UJ 

Ql 



O 


1 o 




u 


o 



Q 


0d 

* 

ad 

O 



OC 











ad 


OC 

* 

a: 

o 



o 


a 




o 


LL 



< 



* 


LU 



o z 


:Z 




H 


H-4 



o 


LU 

* 

UJ 

a: 



LU X 


< 






H- 





z 

# 

z 

OO 



CC 3 


IX 




00 


Z 



LU 


9-4 

# 

*—* 

X H 



— J 


1 X 




LU 


UJ 



X 


H 

# 

H 

o z 



o o 


o 




3 


o 



1 9-4 


3 

* 

3 

< LU 

II 


LU O 


o 




Z 





H 


a 

* 

O 

LU X 


il 

CJ 






i-i 



Csl CO 




cc 

* 

a: 

3 

II o 

II 

3 O 



r 


r 

H 



O (NJ 




GO 

* 

GO 

e> o 

II II z 

II H 3 

O OC 


<SJ sf 

r 

H 

r 

Z 



<M | | 



in 


* 

3 

z a: 

-J LU 

*-H 3 

LU < 


t 1 

r 

1 

r 

O 



* *-l H 



I 

OO 

# 

oo 

< 

J J J H 

H O -J *-* 

OC O 


H CO 

r 

O 

r 

o 






H 


* 








r 


r 








* 









w 










# 

* * 

* 

* * 

* * tt # 

* * * * * 

* * * 

* 

* * * 

* 

* 

* 

* * 

* 

* 

*» * * * 

* 

* 

* 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


CONTINUES TO COLUMN 70 TOTAL OF 12 TIMES 


76 


**************** 


oc 

o 


o 

o 

o 


*** c\ico>frin'Oh-oooo— i cm cn >t in >o r- oo 

* OOOOOOOOtHr-ir-|^r-4r-l— 

* ooooooooooooooooo 

* ooooooooooooooooo 

* UIUJUJUJUJUJUJUJUJUJUJUJUJUJUJUJUJ 

* ccocccccccoccccccccccccccccccccccc 

CC OC CC OC 


o0 

Q 

cc 

< 

o 

LLt 

O 

2 : 

< 

>- 

</> 


h- 



► 


* 

OC OC OC OC OC 

CC UL CC cc cC CC- 

OC QC 

o 



-2 


* 






UJ 



O 


* 






> 



•— « 


* 









*- 


* 






Ui 



< 


* 






2 



u 


* 






i— » 





* 






K 



U. 


* 









►—I 


* 






CD 



h- 

UJ 

* 






2 



2 

o 

* 






i-i 



UJ 

2 

* 






X 



O 

< 

* 






o 



H-l 

1- 

* 






H- 

o 


•» 

00 

* 






< 

• 


o 

HH 

* 






2 

O 


2 

o 

* 







O 


< 

•» 

* 







O 


2 

UJ 

* 






OC 



2 

O 

* 





mm 

o 



O 

2 

* 





o 

1- 



o 

< 

* 





2 

o 



•» 

H 

* 





UI 

UJ 


00 

UJ 

CO 

* 






> 


UJ 

o 

*— i 

* 





•* 


• 

O 

2 

O 

* 





m 

UJ 

</> 

2 

< 

•» 

* 





ff-i 

o 

OC 

< 

K. 

Ui 

* 






2 

UJ 


00 

2 

* 


mm 



ii 

< 

H* 

CO 

►H 

*-» 

* 


z> 





UJ 


O 

H- 

* 





i— » 

CO 

2 

o 

•» 

•> 

* 


-J 



»> 

>-4 



UJ 

UJ 

* 


•* 


<M 


o 

2 

CsJ 

2 

2 

* 


«~t 


rl 

1-4 


>-4 

H 


H 

* 




•» 

w 

il 



H 

h- 

* 


II 


f-l 

H 

O 

UJ 

LL 

•t 

•» 

* 





w 


o 

o 

2 

2 

* 


«~t 


II 

a 

X 

2 


O 

O 

* 


#* 



cn 

h- 

< 

~i 

»—» 

A-# 

* 


mm 


»— » 

H 

o 

H- 

< 

►— 

K 

* 


M 


•» 


2 

l/> 

H— 

< UI 

< 

* 





O 

UJ 


O 

o o o 

* 


3 


>"■« 

< 

_i 

o 

h» 

W 2 

i— * 

* 


M 


w 

UJ 




u. < a. 

* 


— mJ 


1- 

aC 



o 

*-« H- 

i— i 

* 


0 * »• 



2 



r- 

h i/) h 

* 


-J J NC 1 

*■ 

M w 




2 w 2 

* 

— * 

UJ — -1 

H^OfOO 



2 

UJ O UJ 

* 

sr 

l/) CM •* 

•» 

2 

• < 



2 

o •> o 

* 

V CM 

M H Q 

x o uj m uj 



3 

fH UJ 


* 

f-t w 1-4 

•* O h ^ H 

M 

UL OC 



-J 

* 2 

r. 

* 


2 rH 

CM 


CM 


o 

o 


5* O 1-4 
O 2 K 
< < 

oo x 

o 

< o 

t- 

< .. 

o o 


00 
UJ 

— O 

r r 2 
r m r «-• 
r^rK 
^ r i rz 
I rn ro 
h rn r o 
r r 


**************** 


o 

2 

< 


O 

o 


* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 

* * * 


w 


o — 

< in 
UJ W 
CC O 
2 

* * ' 
CC CC 
UJ UJ I 

o o 

UJ UJ 
I- K 
2 2 


o 


O — ’ i 


10 CO 
*-» < 
-J 

a: v 

UJ 2 

o o 

UJ x 
H 2 
2 O 
i-. O 


— . X i-i 
O' H H 

*• •>. *» 

OMQ 

< p* < 
UJ CM UJ 

a: ^ oc 
2 2 
w H ^ 
< 

0 2 0 

< OC < 
UJ O UJ 
a: u. cc 


* o -■* 

cm O cvi 


^tcro 

< O' < 
in • uj 
** <y cc 

UJ 2 
h- • <- 

< O 

X ^ O 

a: < 

OlilU 

ll >-« oc 


* CM 
X H 

in • 

* h- 

in o 

i— i • 

— o 
2 
h- UJ 
< H4 
X w 
OC 

o u. 


co o 

• w 

in 

CL —* 
CM >4’ 


x o 

O < 
«-♦ Ui 
— OC 


< OC 
2 0 0 
OC < H O 
o Ui UJ 2 
U. OC OC LU 


O' r-l 


(NJ 


CO 


O' 

O' 


oooooooooooooooooooo 



77 


r-l 

* 

* 

******** 

* 

* 

* 

* 

* 

* * * 

****** 

* 

* 

* 

* * 

* 

O 

* 


LL UJ 












* 

o 

* 


uj uj a uj x 












* 

OO 

* 


CD 00 O I UJ h 







< 





* 

h- 

* 


00 Q h CO h I 












* 

OJ 

* 


cC OX _J h* LL 







< 





* 

OO 

* 


Mh Joo Q O 







Q 





* 


* 


< H UJ z •* 












* 


* 


Q-ZOC X < Q 00 







LL 





* 


* 


LU *—» LU J~~ 2 2 







O 





* 


* 


H oo < X -30 







QC 





* 


* 


UJ O CL UJ LOOm 







UJ C£ 





* 


* 


00 I I • t— U_ 1— 







>- UJ 





* 


* 


CD O h- O — 00 C 







< CO 





* 


* 


3 Z^-mZh 







-J X 





* 


* 


00 00 LU HXliJ> 







3 





* 


* 


Hi Z J v. UJ UJ UJ 







X 2 





* 


* 


UJ UJ <t 3 CD O 







H- 





* 

++> 

* 


X 2 -J 2 uj 







a I uj 





* 

2 

* 


KOUJO — 200O 







LU Mi x 





* 

Mi 

* 


HH M O O < 0^ 







O H- 





* 

X 



</) h Q h - X < 







QC UJ QC 





* 

< 

* 


h m j < LL O 







UJ CL O O 





* 

M 

* 


OOhhh — mUJZ 







K O LL 2 





* 

00 

* 


uj a; 3 —* t- < 







2 -J Mi 





* 

H 

* 


h <LU a W ►< h 






—J OC 

Mi oo oo 2 


CO 



* 

a 

* 


aXr3«Qw 






LU UJ 

QC UJ Mi 


UJ 



* 


* 


Q K O z _ 






O >- 

LL LL LU QC < 


QC 



* 


* 


2 UJ O ■— Mi O O 






Mi < 

o o >- < ►— 


< 



* 

< 

* 


< X 00 K O 2 2 






-J 

a: <32 


3 



* 

K 

* 


1- mi LU “5 < < 






00 X 

zuizjqo 


Of 



* 

-J 

* 


► Xh-O 






QC X 

O > O U)U 


00 



* 

UJ 

* 


Q LL — H < Q OO 

• 





O 

Mi < Mi X 





* 

Q 

* 


a CD ►OO h 

< 





►— — t 

h JhHlL J 


U- 


K 

* 


* 


O 3 CD — «✓) Z 






O Mi 

< < I OLU0LO 


3 

* 

-J 

* 


2 — w • 3 2 OC UJ 

-J 





LU - 

Mi X Mi Mi O LU 



CL 

* 

UJ 

* 


< CD 3 Mi Mi 

UJ 





> < u. 

> K > X Mi > 

X 


H 

* 

o 

* 


jOhUJOkO 

o 


O rH 


K O 

LU | UJ LL. 3 

< 

3 


3 

* 

M 

* 


1- — 2 X -0 00 Ml 






UJ _J 

Q Mi o O CO X 

_J OO 


O 

* 

M 

* 


D£Chw<H 

> 


O 


oc 

UlUh 

1- 




M 

* 

o 

* 


00 Q O O K UJ LL 

< 


2 


o 

2 o a 

QllQ> JL5I 

-J 


2 

* 

M 

* 


a: zcd hQ -j uj 

QC 


< 3 K 

< UJ 

o£0a:t-<20< 


Mi 

* 


* 


0 3 O 2 00 O 

QC 



Mi 

UhllU 

< <hDI1I< 

3 


o 

* 

m 

* 


KOC^UJ<hOO 

< 


1- 


LU 

oo O QC 

oluooojujo 



* 

3 

* 


O CD UJ > 2 




LL 

> 

Ml UJ 

2 Q- 2 O M, 


Mi 


o 

* 

m i 

* 


UJ CL — Q UJ 

2 


a. o 


O 2 K 

< O < -J OO LL 

2 00 


M 

* 

m 

* 


> a: cl x »-* ll x 

Mi 


o 


I- 

O 2 

H- J H LU UJ O 

Mi 

UJ 


2 

* 

K 

* 


IUDO-QhK 




X 2 

Q M H 

00 to 00 > oc 


QC 


Mi 

* 

00 

* 


OO x OC O 2 

O 


X 1- 

Mi 

2 00 

QC 

O0 



O 

* 

UJ 

* 


2 O UJ LL - < 

UJ 


fr- O O < 2 || 

II II II II II o I- 

— J 


3 

* 

K 

* 


O J X — • • 

QC 


O 2 

a 

UJ 

H* 2 < 


OO 

* 

i— i 

* 


Mi hOHh-- 

O 


2 UJ 


UJ 2 M 


Mi 

I- 



* 

•» 

* 


1- UJ UJ W 

h- 


LU 

-J 


X Mi Mi 

Mi Mi Mi Mi Mi UJ CD CD 



* 

Q 

* 


MIX h h U. h *»> 

oo 


— 1 


oo 

Mi Q M 


1- 


«• 

* 

2 

* 


K H O 0 — 0 




II 

UJ 

h- fH CsJ CO ^ in o 




a 

* 

UJ 

* 


a: 2 UJ K 3 - 

LU 


II 


H 

II — 



II 


UJ 

* 

Mi 

* 


< Mi _| < Z Mi «*» 

QC 



H 


II <<<<<< 




<£ 

* 

w 

* 


Cl o O UJ o — m< 

< 


O oo 

II 

-J h» 

hhhhh*/) 


Z 


Mi 

* 


* 


h O. (/) CO m O ^ 



2 UJ 


O LU J 

J J J J J H 


Mi 


3 

* 

00 

* 


00 z -J K - H- 

00 


UJ H 

3 

— QUJUJUJIUUJUJO- 


X 


Of 

* 

K 

* 


H Mi Z UJ 

K 


Mi 

Mi 

Ml 

HmiQQOOQOZ 


< 


LU 

* 

UJ 

* 


LU O mi aC LU | K 

2 










QC 

* 

00 

* 


oo a. o < x a: mi < 

UJ 











* 


* 


-5 h- < — 

X 










00 

* 

UJ 

* 


uj z oo a 3 oo 

o 










UJ 

* 

2 

* 


ZhqCHH Ml z 

UJ 










2 

* 

Mi 

* 


MiOUJCLUJX— mi 

oo 










Mi 

* 

H- 

* 


h’O^M(^Kh< 











h- 

* 

3 

* 


3 O aC I X 

QC 










3 

* 

O 

* 


0< JUHmoOUJ 

< 










O 

* 

a: 

* 


qC oo oo *-. CC 

LU 










QC 

* 

CD 

* 


CD OC UJ CO QC UJ 

2 










CD 

* 

3 

* 


3OX3M.XC0CD 

M 










3 

* 

00 

* 


LOU.HLOIL WD J 











00 

* 


* 














* 


* 














* 


* 














* 


* 

* 

******** 

* 

* 

* 

* 

* 

* * * 

****** 

* 

* 

* 

* * 

* 


oooooooooooooooooooooooooooooooo 


REAL*8 8ETA(12»23)»STD(8»23) » GAMMA ( 5»22)»DEITA(6»6) SETS002 
REAL*8 TMIN t SMIN t AMIN,SAVE<6l/6*0.00/ SETS003 
INTEGERS NREAD/l/,NPUNCH/2/,NPRINT/3/ SETS004 
REAL T(24» ,D(24I t STC24) ,SD(24I SETS005 



78 


sor^ooa'Ot-H(N4ro^invor^ooO'o»-<rgrO'4“in>or-QOO'Or^(\4rn4-invor , ^cochO— icsj 

OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 
l/)tOCOt/>tOtOtOtOtOtOLOCOtOtOLOLOtOtOtOtOl/>LOLOLOLOCOL/)L/>i/>tOLOL/>LOLOLOL/>l/> 
I— !•— h-H- H-fr— I— I— h-l— H-h-^l— H- h— h-l—l'—H-H'H* 

LULUmLULULULULUUJLULULULULULULULUUULUUJLULULULULULULULULUlJULULULLJUJLUUJUJ 
t/)l/)tO(/)U}(/)OOt/)i/)l/)OOCO(/)l/)l/)tOl^l/)U)l/)t/)^(/)lO(/)lOl/)oOl/)(/)l/)OOcOl/)^V)l/) 


co 


II 

X 

UJ 

o 


00 


o 

m 


Osl 

rH 

o 


5 : 

to 


X 

to 


X 

o 


in 

« 

CO 

-H 

2 o 


X • 

2: 00 II 

>-< •* ►— 


o 

o 


o 

r- 

Q 


o 

o 


or 

iii 


0m 



















h- — 4-4 X 

or 






• 


;2 2 ~ 

rH 



















4-4 IO X 3 

UJ 






Q 


mZmZ 

w 



















•> w tO X 

• 






2 

rH 

X 4-4 5 ; 4-1 

to 



















2 0-4-4 4-1 

2 






< 


F* X F~ X 

h- 



















4—4 »— » O' Z 

4—4 






• 

• 

4-4 f- 4-4 K 

Ql 



















X *-2 4-4 

X 

O 

O 



1 — 

•Ml 


•'4—4 4— 4 

ST 



















K X 4-4 X 

to 


•—4 

O 


X 

4—4 

2 4" •* in ** 

•* 



















► til 51 

w 



CM 


4-4 

0 

4-4 

w in w *o 




















X O H • 

* 

O 

O 



X 0 X 

< ^ < w 

cm 



















UIZ 4 . 4 

Q 

H* 

K 

O 


tO LU |— 

hQhO 

<\J 



















QmZ«I 

2 



H- 


4-4 

• 

4—4 

UJ K UJ h 


< 














A 

•M 



25 0-4—4 4-4 

— < 

O 

O 






00 to CD to 

X 

X 











0 



k- 

H* 


A 

w Q J •* 

in • 

O 

O 

O 


+ 

rH 

1 


LU 

X 











in 



U. 

X 


O 

► oo h - 

• «— «* 



O 





II II II II 

0 

< 










H** 




4—4 

4—4 


to 

O •* 4-4 II 

co 0 

— 

«— » 



1 — • 

+ 

-J 


2 

0 









K 

LU O 



X 

X 


•> 

V)h X 

4-1 r- 

O 

O 






M M M 

4 — i 






a 




to 

4-4 

l- 



to 

to 


•— 

•* to — LU Q a 


r- 

2 


+ 


U 

cO ro co cO 

•» 

0 





2 




UJ X 




4-4 

4-4 


to 

0 ^ »"H 0 

0 - « 

0 

0 

4—4 



05 



— •» 

h- 





LU 




fr- 

to 0 






•> 

0 * -J ► 2 

- ?— 1 

0 

• 

X 


2 < 


+ + + + 

O 

OO 





»— 4 




4— » 

4-4 

O 



+ 

+ 


-J 

_J w H-. 

11 • 

#— c 

rH 

CO 


4—4 

X 

CO 



•* 









4* 




-J 




w 

w 2 

to 0 

• 

• 

0 


X 



J J j j 


< 





II 





I 

*** 


•» 

-0 

-0 h- 


~ K 

UJ 

a 

a 

hUih 

w 

+ 

LU LU LU LU 

Z> 

H 


0 









4~ 


r*H 

w 


LL 

2 

2 Z> 05 

X • 

LU 

LU 

0 h- 

4-4 

w 


OOOO 

M 

UJ 



O H 

r— 




II 

■ 0 * * 

• 



h" O 

4-4 

4—4 

m a a. • 

0 2 

0 

• 

0 

< 


• 

-J 

H H M H 


CD 

O 

It 


I 

h 





4—4 



H 



X 

O 

a 2 ~ 

X v-i 

2 

2 

2 a 

H 


LU 

• 4 • 4 

cc 




II 


to 




»—4 

w 

-J 



II 

II 

to 

“5 

-0 O W 

4-4 ?! 

4—1 

4-4 

4—4 

4—1 


O O 

rH (M CO 4“ 

LU 

2 

II 

fr- 


II 

UJ 

O 

H 


X 

• 


-0 



4—4 

to 

OO h- 

2 1- 

X 

X 

X 

O 

h~ 

2 

4—4 

W r Sim* w 

O 

O 


LU 

LU 


H 




0 H-t 

-J 

0 


M. 

M. 



LU < 

4-4 w 

F- 

to 

h- 

2 

LL 

0 

w 

< < <t < 

LU 

X 

— 1 

4—4 

> X 

4-4 

II 

1) 

II 

in 


w 


00 

-0 

-O 

II 

-J 

— 1 — 1 K X 

X W 

w- 

w 

w 

< 

4-4 


t/)K hhh 

h- 

X 

UJ 

X 

< OC 






II 


II 


W 

•4* 


-J 

j j m ot: 





O X 


H- 

J J J J 

2 

O 

O to 

to < 

X> X 

rH 

CNJ 

0 


LL 


0 1-0 

to 

< 

< < a: 0 

• LL 

LL 

X 

u. 


to 

LL 

a 

LU UJ UJ UJ 

4-4 

O 

4—4 

1—4 

4-4 


4—4 

4—4 

4-4 

4—4 

0 

-1 

4—4 

X O to to 

4-4 

O 

0 0 X LL 

4-4 

4-4 

4—4 

4-4 

f— 

4—4 

4-4 

2 0 Q Q 0 


00 


O' 


* 



79 


ro^-LOsor-oocrOF-ic\iro>^tnvOh-cx)crOrHC\iro4'msDr w coO N 0»-HfMrr>'4'U^sOh-coO N 

4 *' 4 - 4 -Ni“>t^- 4 *inir\unininirkLnLnLnLr»s 0 s 0 N 0 v 0 ' 0 'Ov 0 ^ 0 s 0 sor w r s -r^r^-h~N-f^r fc -h-r^ 

ooooooooooooooooooooooooooooooooooooo 

l/)(/)t/)(/)C/)CO^l/)t/)^</)O 0 (/)U)(/)t/)O 0 </)(/)(/)(/)(/)(/)lO(/)l/)(/)t/)l/)^</)t/)^)^l/>^^ 
►-I— K-H-KH-h-HH-H-h-Kh-h-KH'l— h-h-h-f-KH-l— H-h-Kh-H-Kh-H-KH-H-t— >- 
llJOJUJLUUJUJUJUJUJliJUJUJUJUJUJUJUJUJUJUJUJUJUJUJLUUJliJUJUJLUUJUJUJUJUJUJLU 


»-» 3 l 


p* sO 


a: 

a 


oo 


O 


<M 

mm 

< 

h- 

UJ 

CO 


o 

CD 


2 

M 

X 


o 

in 

o 


o 

CD 


2 

v- i 


CD >-4 


2 
v— » 
£ 


in 

•kkkT 

< 


I 

z 


X 2 

*- •-« 

•-f x 

— t— 

i 2 : ~ 


* z 

r H -4 

« z 

X Z h 


w 2 HH 


J X h h 


<N (\J 

1 CO 


ro co <t 

to 


CO CO 


a + cd _i 

2 UJ 


2 : *- 


a. 11 u ii 11 


oc 

o 


2 

♦— 1 00 

X 

CO 

0 

•> H- 

CM 

w a 

< CD 
X 

X -* 

< -> 

o 

b~ 2 
CD ►“* 
• X 
— » CO 
2 

►-* •» 

2 ^ 

H t/) ^ 

CO ►-» < 

- X 

0 4-2 


CO 

W, 


in uj < 

II hh 1“ + + 

+ 

K w < 

V 

< 


W CD h- 

*• + + O 


< O 

0 

I-" 


< • co uj 

iH •* H r— 4 in 


O 2 • 

1“ 

UJ 


H— •— 1 CD 

*-4 _ 1 M (Nl M (VI (\| M 

UJ 

CD 2 I- 

00 

OQ 


UJ CD • 

UUw»-4wl-»H-|W 0 

a 

< CD 



mm — * ■ 

CD * O h- 

+ Q<W<S<W<Nj-h 

*— » 

— * 0 • co 

II 

II 

UJ UJ 

^ * 1— -J 

• H* O 1“ O O H* H 


*— t ^ CM 


UJ 

> > 

^ (O • 

ZwLJJh-UJhhliJ O 

11 

• w m 

— 

— > 

< < 

• O — 

i-hCOCQCOCDOOCOCOOCD 


UJ • 0 


2 0 + 1 - 

• 2 


+ 

+ •-« 






< 

-J 

-A 0 


1—4 

2 







O K 





O 


—k 

< 

u 0 






-J 

+ + 


• 

UJ 

O »-H 

+ 



— * 


«N> 

A9 

mm. 


CD CO 



2 




H 

• 

UJ 0 

+ 

-j 

— 1 1 



9k 



Mb 

O 

t 9» 


n 

A*. 

.J 

-J 

-J 

-J 

— J 

-J 

UJ 



M 





— k 

0 



UJ 

UJ 



H 

-3 -0 


H 

*— 1 

a m 

— J 


rH 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

— H 

H 


9> 


— • 

UJ 

+ 

r-l 

1—4 


-J 

0 

0 H* 


UJ 


kw> W» 


t 

«*> 

UJ w 

UJ 


0 

00000 

O r-l 



«— 


i-H 

H 


• 

9k 

•—4 

UJ 

♦— 4 

h-» UL 

H 

> 

Ii 

a 


O CO 

• < 

0 

-J O 

*-4 

*-4 

H-4 

1—4 

►— 1 

»-4 

• • 

+ 


CO 


• 

< 


O CO 

• 

Q 

► 

9- M 


< 




UJ 

w 

— 1 K 

*-4 

UJ 

UJ 

9> 

► 

9* 

9* 

*» 

9> 

O K 



i-4 


O O 

O 0£ UJ 

k— " 

a 

*—4 

in 

0 X 

+ 

00 

T> 

11 11 

O' 

• 

< UJ -J 


0 

• 

H 

C\J 

ro 


in 

O 

UJ -J 

2 




uj in 

»— » 

< 

• 

< 

UJ 



^ 00 


*— * 




X 

1“ 

O UJ 

II 

HH 

2 

w 

w 

w 

«w> 


w 

• • 

•— 1 

>4* in 


O 

• 

O 

2 

X 

1“ 

• 

II 

< 

< *"H 

x 


r-4 

«k r*% 

0 


-J 

>■* O 


w 

>— 1 

<£ 

< 

< 

< 

< 

< 




O 


* O 

2 

w 

k— - 

-J 



»“ 

H* 


11 

H 

-9 -9 



UJ 

W- W 

-J 

00 

w 

K* 

K 

1“ 

1“ 


H 


II 

11 11 

1“ 

II 

^ I- < 

■«— ' 


UJ 


_J 

-j 

-I II 

ii 



w 



0 


UJ 

H- 


_J 

«J 

-J 

-J 

-j 

-J 







O 



0 


UJ 

UJ 

UJ 


CO 0 <— 0 0 LL 

mm 

IL CL 

O CL 

CL 

UJ 

UJ 

UJ 

UJ 

cu 

UJ 

UL CL 

J ^(MOZILO 


LL 

LL 


u. 

O 

O 

O -J 

x 

►H 

Q IO 00 CD 

M 



1— 1 

2 

HH 

Q 

0 

O 

0 

Q 

O 

(—4 »-H 

1—4 

M M 

CD 


»— » CD 00 

*— t 

»-4 


w 

*— 1 









X 




















X 










CM 



CO 










st 






O 










r-l 


H 



•-4 










9-4 






CM 






O 



80 


QrHOsjrn^m>or^ooo^OrHf\iro>4-mNOf s -ooo'Or-4<\jrOstLnvOf s “aoa s o*-<(>JfO'^ , ^'0 

COQOOOOOCOOOOOQOOOQOO'O'O'O'O'O'O'O'O'O'OOOOOOOOOOHH^HHH^ 

OOOOOOOOOOOOOOOOOOOO.H.H*HfHrHrH.HrH.H,H.H.H*H.HrHrH,H 

t/)U)U0COU)COWO1tOCOOOt/M/)</)t/)t/)t/)V/)(/)(/)VM/)i/)t/)^CO(/)(/)U)l/)OOCO(/)O0(/)t0t/) 

LULULUllJUJLULUUJUJLUUJLUUJLJJLUUJUJUJUJLULUlUUJUJUJLUUJLUUJUJUJLUUJUJUJLULU 

00O0o0o0^o0o0U)^^0OU)t0^C0COO0c0^tOl/)i/)(^</)O0tO00^U0</)U)O000(/)(/)00(/) 


UJ 

a: 

<£ 


50 

o 

TO 

trv 

O 

GO 

H- 

O 

0 

CD 

z 

— 


00 

K 

z 

UJ O 
*-< a£ 
o <i 

4-* O 

IL Z 

a. < 

UJ K— 

o oo 

u 

UJ 

ui o 
*— 0 C 

< H- 

O 

1-4 UJ 
O X 
Z H- 

< 

O U. 

O 

00 

00 
LU Z 
Z X O 


X 

00 


00 


X 

00 


CM 

< 

X 

X 


a 

*H 

o* 

CO 

< 

h- 


X H 
00 OL < 
H-» O X 

1-4 

I 00 X 
Z O 

jOfli 

i— « CL 

II h ft* 


•-« — + 

X Z 

— 00 M H 

Z »h X >-♦ 

»-H — ^ 00 ^ 

X Z .H 4-. < 
oo W ► X 
•-* X 4* #H X 

•> 00 < 



H 





— 

— 

CM 




O 


< UJ 




< 

< 


rH 4-4 

rH 

♦ CD 



— CM 





z 

z 

w » 




o 


O ro o 




— 4-4 



4—4 #• 

4-4 







LL 


♦—4 — 

4—4 

< z 




• 


• CM • 



H 

-» > 

LU 


w (M 


CM 



w o 



*—» 


X z 

X 

Z X 4-4 




o 


H- f— 




UJ UJ 

a: 


< «H 

< 

4-4 



UJ H- 



X 


CO 4-4 

00 

»-» X X 






CD O CD 



+ 

a a 

< 


X w 

X 

w O 



> 



00 


4-» X 

1-4 

X < 00 




II 


• |- • 




t-4 



X o x 

o a 



< O 



4-4 


► </) 

•• 

00 CD W 


o 




— — . 



z 

w Q 

o 


< fr— < H- • 



oo CD 





rH «-« 

CM 

H-4 •» 






z o z 



1-4 

00 CC 

z 


CD 00 CD 00 rH 


UJ 




+ 


— * r* 

W 

•‘'SCO 


II 


CO 


►h CD *-< 



X 

H- < 

< 






to 

U — * 




H< H < 

CM — 






X X 


HWQ.Q 



U II 

II 

u ii 


H» 1 

K 


H- 

Z 


X - 

X 

w o < 


A 


+ 


00 — 00 



4—4 

z z 

o 







o — oo 


LU 

•— i 

+ 

X O X 

o o 


m «o 



4-4 O 4-4 


+ 


< 

H* 


— ' •— k 


— — 


It 

► J UJ 


4-4 

X 


< K 

< 

K • UJ 



•* 

«J 





II 

— f— 



-J 

— J 

-j _t 



HUJH 


X OO 

ZOWOWHCQ 


+ 

*— 4 

UJ 


>* C* (\J 

* 

•J 


f-4 OO 


• 

UJ UJ 

UJ 

UJ UJ 

rH 

A 

Q <-» 


OO 

4—4 

4—1 







Q 


w UJ — 


UJ 

— 

• 

X to Q Q 

o o o 



II m • 


4-4 


X 

II II 

II 

II II II 



II 

4—4 


< • < 


o _j ar lu 

O 2 h h 

4H 

4-4 4-4 

+ 

UJ 

•* Of M 


II 

00 




o UJ 


•* o 

X -J X 


4-1 

LU 

UJ x 

ot O 

•• •» 

•* 

•* »• 


Q 

*o *) UJ 


II 


4-4 


— k 

— — — . 

in o 

-o 

-> m 

x uj x 



Q 

• ►— 

U- 

4-4 

rH CM 

CO 

sf in 

X 

4-4 

V • 

» 


h- 


0—4 CM 

ro 

4 

5 

6 


•—4 


w 


< o c 

0 

l 

l 

II 


bd 


K 

<w w 


w w 

X 

w 

O <t 4-4 


UJ 

u. 

II 

— ^ w 

«— i * 

w w ^ 

o 

w 

O' 

< 

Q O w O 



w 

b£ •• 

z 

< 

< < 

< 

< < 


oo CO t— ^ 

b£ 

> 

♦—4 


UJ LU 

UJ 

UJ LU UJ 

h OO H 


K 

W W 

II II II 

_J 

oo 

W LU 

UJ 

4-4 

H H 


I— 1 - 

II 

K 

—1 

a: 

< 

X 

UJ 

> > 

> 

> > > 


h- 






UJ I-. 

1 - 

X 

> 

-J -J 

«J 

-J — 1 


a 

□ UIIL 

< O 0 00 

to 

< < 

< 

< < < a 

Q. 

O UJ O 

LU LL LL 

ylHovjo 

a 

IL O 

<f UJ UJ UJ 

UJ 

UJ UJ 

b^ 

z 

O O *-H 

X 

t-M 

*— 4 

1—4 

00 00 

oo 

oo i/> 00 CD Z 

Q 

Q CD H w H 

b£ ~ 

»-4 

z 

»-4 2 

1—000 

a 

o a 

bC 


O H 

ro cn 


a) O' m cm 

rH »H CM CM CM 


u u u 



81 


r^co^or-*c\jm>*mor~aoo'0»H<\irOvfrLnsor^coo>o«~i<\J^>J , inoh“ 

rHr-*«-»csJc\jcMc\j(Nj(\jOsirsJro(\jrornrOfOrornrOfOrorO‘>i , 4*^^'>t>t'4*4“ 

h^hhhhhhhh^hhhhhhhhhhhhhhhhhhhh 

COcOtOLOtOtOtOtOtOCOCOCOCOCOCOtOcOCOCOCOCOtOCOCOoOCOtOtOCOtOtO 

UJUJUUlilLUUJLiJliJUJLUUJLLJLJJLiJLUUJUUUJUJUJUJLUUJLUUJLUUJUJLUUJaJ 

tOtOtOCOoOOOCOtOtOtOCOCOCOtOCOtOcOCOCOtOtOCOtOCOtOCOCOcOtOtOCO 






z 



UJ 



to 

• 


o 

h 


X 

LU 


o 

CO 




• 

to 

ID 

Z 

1— 1 

K 

HH 


UU 

X 

H 

-J 

CO 


Ol 


X 

x 


3 

o 

z 

X 

o 

M 



X 

z 

LU 

to 

%-< 

X 

p» 

X 

h- 

z 




o 


o 


u o o o u o 


3* VELOCITIES MUST MONOTONICALLY INCREASE. SETS148 

SETSI49 

END SETS150 



82 


f-l 

45* 45* 

4ft* # 

* 

* * Hr 

4ft- 

* * 4ft- # 

* 

* 

4ft-4ft-45-4ft*45-45*4fr4fr4l*4ft*4ft“45- 

* 

* 

eg co 

^ in 0 

o 

* 










* 

0 0 

000 

o 

* 



eg 





00 00 


* 

O 0 

000 

2 









LU UJ 


* 

2 2 

2 2 2 

»-* 

* 



* 





QC DC 


* 

M HH 

H M H 

0 

* 








< < UO 


* 

O O 

OOO 

3 

* 



2 





3 3 UJ 


* 

3 3 

3 3 3 

1/) 

* 








O' 0 z 


* 

00 l/) 

00 00 CO 


* 



II 





CO </) w 






* 








—1 


* 




4ft* 

H 


MH 





LL LL 






* 

*“• -J 







OOO 






■* 

LL ULi 


O 





2 


* 




45* 

O 


ft- 





X X ft- 


41* 




* 

H- O 







3 3 






* 

00 2 


eg 





CO 00 LU 


* 




* 

OJ 







X 


* 




45* 

CD OL 


ll 





a a k- 00 


* 




* 

<C 







»—* »— ♦ LU 


* 




* 

LL UU 


•— « 





00 UO LL QC 


* 




* 

O 2 







UJ UJ 0 < 


* 





•— 1 


X 





OC cC 3 


* 




* 

00 -J 


0 





z 0 


* 




45 * 

2 


ac 





II ll O 00 


* 




* 

O UJ 


LL • 





H-4 






* 

>-< -J 


— » X 





H-4 4 1— LL 


* 




* 

t— CL 


2 00 





•* * O O *-i •-< c\ieg 


•fr 




* 

< 2 


W UJ J 





CO SO UJ < D h < (D h 


* 




* 

2 *-h 

*-4 

X 2 _J 





X 2 2 


* 




* 

*-* co 

w 

^ < 





< < c£ 3 LL UL UJ LL LL UJ 


*■ 




* 

CD 

XV- j 


oc 



hhUJtOOQw OO^h 






* 

X < 3 : 


a: 


0 



UJ LU H* O O 


45* 




* 

O 

V 

X LL 0 


OC 1- 



C0CD2J22W 2 2 h-4 






4ft* 

O 2! 


0 LL 


0 0 



«<OOil O O u. 


* 




* 

O 

X 

V 


H LU 



♦» •* DhhILHmmILcM 


* 


0 m 


* 

LL 


00 X 


O > 



H Csl Ui O h K UJ CD h h OJ CQ 


45* 


eg 


* 

a ci v 

— 2 


LU 



GO CO X m < < a < < 0 


* 


eg 


* 

OJ 


rH O H 


> UJ 



(1 ||*-00mh(JLLmmOLL 


* 


•» 


* 

00 00 

—* 

»— 4 CD 


-J 



— ^ LU > > 0>> O 


* 


in 


45* 

me 

t-4 

1-4 K 


LU CO — » 



m m u. a: uj uj z lu ui 2 


■H* 




* 

hh qq 

w 

~ < + 

> 

— 1 < H 



► * 0 QQ030Q0-J 


* 


< 


4ft* 

Q£ 

X 

X 2 


COh 1 - 



eg in Q »-h < < 


•K* 


X 


45* 

UJ UJ 


H-4 r-t O 

< 0C2H 



^^LUUJOOKOQOHU 


41* 


X 


* 

00 oc 

OC qc CD < 

2 w < ~ \ 



<<320CQC<OQCQC<O 


4f 


< 


* 

< C> Q X 

< 

QC > Z 


•• 

hh jm<< ja:<< 


45* 


0 


45* 

< 

LL 

LL O II 


< X *~* 


00 

LULU <CDOQUUQ. QQUJa 


45* 


•» 


* 

X 


O 

X 

> H- 


<coco>x22a:»-422a:*-« 


45* 




■* 

00 O 


> 


Z N X 



O <<a:o<<a;o 


* 


CO 


4* 

UJ »~4 

X 

X LU 

LL 

KlilH 


0 

- ^XOKKOUihhOUi 


45* 


eg 


* 

H* X 


OO J O 

z 0 00 


LU 

H C\J OOOOOQCOOOOUOC 


45* 


•• 


* 

< 2 H CsJ O UJ 


OJ 2 O 


a: 

< < 11 11 


45* 

KJ 

00 

X 

# 

-J 

GO 

CD X O 

X 

Q LU LU O 


0 

II II II II II II II II II II 


45* 

1 

w 


4ft* 

3 0 . 


ft- a 

h~ 

2 CL 2 UJ 


H* 



45“ 

O 

0 

> 

45 * 

O X 

+ 

+ X 

O 

LU UJ O 2 


00 



45“ 

•» 

h 

»> 

* 

— 1 


O0 

2 

CL O h O 



•* •» •» M-4MHHHWHM 


45* 

X 

OO 

2 : 

4ft- 

< 2 

f-H 

eg lu oc 

UJ 

UJ Z 00 H 




45* 

1 

•* 

w 

45* 

O O < 

< 00 < 

3 

O H z 00 


Q£ 



45* 

< 



* 

CL 


►—» aj 


UJ 2 




45* 


CO 

2 

4ft* 

2 3 

II 

II 0C 2 

II 

II II X UJ 



h-KKHOQQOOOOO 


45* 

OC 

eg 

M 

■fr 

1 — < 


Q- >-< 


X 


< 

UJUJUiLU|-hl-Hhl" hh 


45* 

00 UJ 


O 

* 

O >- 

> 

> X J 2 

> X Q hH 


H 

CDCDCDQDOOOOOOOOI/JOOOOOO 


4fr 

45* X 

— * C\J 


45* 

3 


O 


O 


< 



4t 

— J 3 H H 

00 

4t 

00 LL 


U UJ 


UJ 


0 



45* 

< 2 

W> w 


* 

O 


~J 


CO UJ 





45* 

UJ •> 

> < 

LU 

4 * 

UJ 


00 a 


CD 


H 



45* 

CC X 

4— • 

2 

45- 

2 00 


LU X 


ft— 


3 



45* 

3 

— UJ 

HH 

4ft- 

*-* UJ 


t-H HH 


00 1- 


a 



45* 

►— 2 

r-H CD | 

5- 

45 * 

h h- 


a: 00 


3 OO 


H 



45* 


w 

2) 

* 

3 < 


UJ 


X 3 


3 



45* 

O CO 

X 2 2 

O 

45* 

0 2 : 


CO < 


X 


O 



45- 

*-* 45* 

O 

q: 

* 

a: >-» 




< 





45- 

-J -J 

-J X II 

CD 

45* 

CD K 


LU O 


ft- 0 


UJ 



Hr 

a. < 

< X 

3 

45* 

3 00 


X 2 


UJ H 


X 



45* 

X UJ 

UJ 0 3 

00 

45* 

00 UJ 


h <£ 


CD 00 


h- 



45- 

QC 

OC O 3 


4ft* 










45“ 




* 










45- 




4t 










45- 




45- * 

* Hr 

# 

4ft- * 4ft- 

4ft- 

* * * * 

* 

* 


* 

45“ 




ooooooooooooooooooooooooooooooo 




83 


r-oocr'Of-<r\imsi-u'\'Or^aoo'Or-«cvj{'ONfrmvor'-aoo'0— ic'gcn-j-u'i'Or^cocj'O.-tovioo 

o o O h h HHi-tHHHHNNNMNNcgNNNi'nfomfnooforommm't't't't 
ooooooooooooooooooooooooooooooooooooo 
Z2ZZZZZ2ZZZZZZ2ZZZZZZZZ2ZZZZZZZZZZZZZ 

HHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHH 

ooooooooooooooooooooooooooooooooooooo 

^Ootnooooc^oot^oocoooi^oooovjcoooooi/jcocooocototocoooootocoooco^^tocoto 



o 


84 


<}-ir\sor~ooa-Of-icMm^ir\'Or'-ooo'Or-»cMc^'i-inyor'-ooc r 'O , ~ic'jon'frin'0r^<x>O'O 

•4-'4-'4--4-'t'j-irvir\trvu n vir\ir\tnir\intnvOs0vOsO'£ivO'O'OsO'JDr-r~-r-(^f^f^f w f v -f v -f^® 

ooooooooooooooooooooooooooooooooooooo 

2Z22Z2ZZZZZZZ2ZZZZZ22ZZZZZZ2Z2ZZZZZ2Z 

ooooooooooooooooooooooooooooooooooooo 

</>00</)COOOC/300uOl/)COL/>OOOOOOCOl/>l/>OOt/)COLD</)COt/)COCOCOOOOOl/>l/>COOOOOOOOO</> 


in 

r-4 

a 





«•*» 

O — 




2 . 







X 

o o 

r- 










L/> 

2 

rH 



> 







* 

— LU 




8 







X 

— Q 

o 










co 

O w 

1— 



H 







+ 

Q h- 




1 





«*** 


2 

• a: 

o 



2 





2 


O 

O Of 

o 








O 

CO 

2 

• 00 




X 2 





2 

rH 

UJ 

LU Q 

A 



* ~ 





LU 


Q 

-J * 

o 



~ X 





O 

o 


• — ■». 

in 



-0 * 





«—• 

H 

* 

o 2 

i 



*» «— • 





* s 



2 o 

o ~ 

CO 


eg -> 



A 


o 

O 


LU 2 

• -5 

CNJ 


w »» 





?: z 

o 

O 

Q LU 

rH •* 



< tn 



#• 


Q LU 


* 

~ Q 

o eg 

o 


H w 


«0« 

in 


2 O 


o 

• w 

i~ ~ 

K 


LU < 


2 



LU 

o 

w 

QC K 

-i ~ 



CO K 


»— «• 

< 


Q * 

o 

V 

O OC 

• < 

O 

o 

+ UJ 


X 

H* 

A 


e 

CD — . 

• a 

— * K 

O CNJ 

—• CD 

"0 


Ul 

“5 

* O 

O 

•w CD 

00 

— * UJ 



T> 

•* 

# 

CO 

•* 


• 

'■w' -w 

o a 

*0 CD 

— * 

o 

- 1 

CO 



CO 



O 



h- 


K 

K O — 

eg * 



o 

— <• 

H* 

rH 

w <r^ «m> 

V. 



w o 

o 

v. a 


x u> o o 

a: 

a: ♦ ^ 

O CM O 


I s - rH 



O < 2 2 




o < o 

a 

• 

X 

X 

• 

Q O 

a a o 2 

l-« *1 O 


o 


O < 2 

O h* w 

o o 



O K • 

• 

< o oo 00 

* 

• 0 

oo 

oo • D 

• w 

t 


• 

l 

Oh- — 

• LU X X 

a o 



• LU O 

eg 




-o 

o o 

o 

O LU 2 

o < 

rH 


rH 



LU > 

O CO 

I s - • • 



h- O CD 

i 

11 11 

II 

It 

•* 



_j 

K- 




•o 


CD 

11 II 

►Oh 



* II 

o 




CO 

n ii 

ii 

ii • it 

II LU 

II 


II 


-o 

— II 

II II 

in 



^ II II 

w 

-•>« 



»— 4 



2 

CD 




w 

-o 

II 


M II 




w 

-5 

"0 

-o 

w 


«*»» 

— O -* 


<**> 


«•» 

t 

• 

a «■* 

— *0 ^ 

II 



II — > — -o 


•* eg •* 

•* 

<* <t 

oog 

-o 

OZOnOOOO 

-o 

00 


l- Q> “5 

-o -0 ► 



in 

-o *> *• 

N» 

CO rH CO 

H CV| h 

•* rH 

► 

•» UJ ^ 

rH •* CD 

•» 

rH 

r. 

_J 

UJ 

•* •* 

•- •* H N 

HH -0 -0 


(VI 

HH •» •* »H 


>— 4 HH 

*0 

-o 

UJ 

rH eg 

HcgQfn 

ro < <J- 


4* 

• 

• 

in 4* 

vO CD rH rH 

•» 

2 


m 00 H 

< 

~ 0 — 



CD 

X X o 

X 

X X 

O X Q 

X 

o 

X 

-o 



W- W W w 

O' ^ 00 


O rH w «“* 


< H < 

< < 

w 

H- 

w 

w W 

K ^ w 


H- 




< < 

< < < <C *H w ^ 

II 

H 

CM < < C 

11 

H- 1— H H 


Q O 

a o a 

O 

a 


a 



y— h- 

1— H" H- L— 

o a 



f— H* t— 


UJ O LU 

LU 

UJ 

LU 

1— h- O K- 

I-ILH-OK-LL 

*- 

Ol-U. 

u. 

LU UJ 

UJLUUJLU0HK200UJUJLU 

CD 

CD O CD 

CD 

CO hh 

CO CO o 

oo 

(/)huU5(/Im 

00 O 00 



CD CD 

CQCDCOCOOOOOO-OCDQCDCDCD 


rH 

eg 




CO 


in gj 



e- oo 




ce 



O rH 


»H 

«H 




rH 

rH 

rH rH 



rH 

rH 




iH 



<M eg 


o< 


85 


*H(\jm» 4 “mvOf ,fc 'Ooa'O*-M\jro> 4 *Lnv 0 r t -oo 
oooooocooooccocococro'a'O'a'^O'a'O' 
o o o oo 0000000000000 

2Z2ZZ2ZZZZZZZZZ22Z 

hhhhnhhhhhhmnhhhhh 

000000000000000000 

OOOOOOOOOOOOOOOOOO 

t/>t/) 00 </) 00 </) 0000 (/)</)COC/>UOl/)</)t/>l/)l/> 



0* 

O 



2 


0 

< 


•» 

• 


in 

— 





< 

a 


►— 

X 


UJ 

O 


CD 

w»- 

0 

1 

X 

CM 

— * 

• 


O K 

O 

•* 

0 

H- 

CM 

* 




O 

< 

O 

O 

h 

*» 


uj r- 

tmm 

CD 


O 


< 

O 

N* 

H O 

0m • 

0^ 

uj r- 

O O 

0+* 

CD Q 

*» • 

0 

W 0 

O h- 

► 

• H 

— CD 

rH C* 

< * 

w 

O II 

h- — » 

< 

• 

UJ 

1 — 

Mb Mb 

CD O 

UJ 

— ■ O 

*> 

GO 

f-» * 

+ in 

1 

1 00 


2: 





— < 

O O < 



CM 

O K 

’ h- 




♦. UJ 

^ X UJ 



O 

CO CO 

w • CD 

0 


1 — 

w O 

< K 

0 



< 1 h- 

1 — -J 

• 

0 

O 

K O 

LU • — . 

0 

0 

CD 

UJ — • 

CO — — * 


00* CD O H w O — » 

It ► O H 

in HH II CNJ H II h- | 

— II • w O 

Oil LU4* — < — — <0 

•* -flO Z OKOr^OH^LU 

(\J M o N ♦ Z O' ► UJ •* Ovl ** UJ • o 
r-4 ► -lOZ CO CD f- NCDHZZ 

cm ►-< o »-< o^s-^o^^jwq: 

< OJ ^ h ^ II II h < *-• < h < • I- O 

1-0 J-I-K^OZKQ 

iuohOita:zouuu.ujouju.^oiiJZ 
cqq woh-?oO(Dn(docDn oo:uj 

X 

cnj ro >fr in or- 00 

CNJ CNJ CNJ C\J CM CNJ <NJ 



86 


H 

* * 

* * * 

* * * * 

******* 

* * * * * 

(M«l4 , in'ON®0'0^iMf)'t 

o 

* 



Q — 

* 

O O O O O O O O r~4 1-4 r~4 r-4 

o 

* 



Z — Z 

* 

000(1)000000000 

z 

* 



< z «-« 

* 

ZZZ2Z2ZZ2ZZZZ 

*— 4 

* 



4-* X 

* 

hhmhmmhhhmhhh 

o 

* 


• 

•* X oO 

* 

ooooooooooooo 

3 

* 


O UJ 

K K h-< 

* 

-5-)~?-5-5-3-3-?-5“3-3-5“3 

u 

* 


H O 

h- 3 ~ w 

* 

OOOOOOOOOOOOO 


* 


< UJ 

< o 

* 



* 


a x 

X z X 

* 



* 

H- 

UJ < 

hOMUJ 

* 



* 

•“« _l 

H UJ U 

UJ o Q 

* 



* 

LL UJ 

-J CD 4-4 

LU ^ Z 

* 



* 

O 

CL 000(/)h 

* 



* 

h“ O 

Q-h- Z 

O 

* 



* 

oo x 

<</)*-! 

Z J X 00 

* 




* 

UJ 



3 


a 

O K 



00 

* 














* 

GO 

OL 


UJ X 

o 

4-1 UJ 

£j£ 4—1 



UJ 

* 














* 


< 


CC 

H 

h- a: 

u. 



CC 

* 













— * 

* 

LL 

UJ 


< to 


< < 

O 



< 

* 













z 

* 

O Z 


z 

O 

-j 

00 z 



3 

* 













•-» 

* 


1—4 


00 O 

LU 

3 00 

00 < 



O 

* 













X 

* 

00 



h- 4-4 

> 

o i- c£ 



00 

* 













00 

* 

z 



00 H* 

UJ 

-J oo 

mm 




* 













4—4 

* 

O UJ 


LU < 

X 

< UJ •£. Z • 



LL 

* 













#> 

* 

M 

-J 


I- — J 


O K 

3 •-» -> 



O 

* 













z 

* 

1- 

CL 


3 

oo 


X x CM 




* 













4-4 

* 

<t 

X — > 

— 

-J o 

4-4 

UJ -J 

4-» CO 



X 

* 













X 

* 

z 

4—4 •—» 

Z < J 


XJZ - 1 



3 

* 













to 

* 

4-4 

OO 

w> 

QC < 


H- < 

1-4 


H 

CM 00 

* 













•* 

* 

GO 

X 

X 

UJ o 

CM 


X z z 


CD 

CO 

* 













z 

* 

X 

<?; 


> 


Q - 

4-4 *— * 



-J 

* 













hh 

* 

o 

V 

V 

UJ UJ 

1 

UJ H UJ O 


It 

11 < 

* 













X 

* 

o 

z 


oo X 


*- 

x -o x 



3 

* 













b- 

* 


a x 

X 

h- 

Z 

< II 

K O 


mm 

— a 

* 




0 m 









1—4 

* 

LL 





z 

in 


4-4 

4-4 4-4 

* 




CM 









•» 

* 

O 

O V 

V 

• X 

II 


X 


•* 

» 00 

* 




CM 









z 

* 


UJ 


Z O 


x o 

# o o 


CM 

UJ 

* 




•» 









4-4 

* 

00 

OO 

mm 

4—4 

X 

l-H O 

OcCUJ 


W 

— CC 

* 




in 









X 

* 

UJ 

< 4-« 

4—4 

X X 

1— 

-J — J 

UJ LL Z 


< 

< 

* 













h- 

* 

4-4 

CO ^ 

■— «» 

3 

o 

LU 

K O 


X 

X o 

* 




< 









•» 

* 

CC 

X 

x O 

z 

LL 

<1 00 4-1 


X 

X UJ 

* 




X 









X 

* 

UJ 

UJ 


z z 

UJ 

> 4-4 

-J 00 00 


< 

< z 

* 




X 









UJ 

* 

00 

CC CC 

oc 

uj o 

-J 

-J 

3 CC Z 


O O 4H 

* 




< 









Q 

* 


< o 

O -J 


U 

O LU 

•• 


GO 

* 




o 









z 

* 

< 

LL 

LL 

oo 


3 • 

-J UJ X 

00 

•* 

X 

* 













4-4 

* 


X 


to _J 

X 

U- H 

<1X4-4 

< (M O 

* 













► 

* 

oo 

o 


< UJ 

LU 

oo oo o 1— o 


< 

< o 

* 




CO 









X 

* 

UJ 

M X 

X 

X o 

Q 

00 UJ 


o 



* 




CM 









m 

* 

H* 

X 


o 

Z 

UJ H 00 O LU 

UJ 

II 

II II 

* 

— 












X 

* 

< 

3: r-i 

CM 

X X 

4-1 

O 

h Z CD 

CC 



* 

fsj 



00 









o 

* 

-J 

CO 

CD 



3 X 

< 

o 


mm mm 

* 

1 












o 

* 

3 

0* 


CC LL 

CC 

00 U 

tO h— 

K 

4-4 

4-4 4—4 

* 

o 



o 









-J 

* 

o 

X + 

+ 

o o 

o 

4-4 

-J — * OO 

00 

»> 

•* #. 

* 

•» 



h* 










* 

— J 



K 

h- 

00 X 

UJ Z 3 


H 

co in 

* 

X 



00 









z 

* 

< 

Z r-4 

CM 

O Qd 

U 

< 3 

O *-» X 

LU 

»>— > 

— ' ’ w 

* 

1 



m 









««• 

* 

O 

o < 

< UJ UJ 

UJ 

X 

O X 

qc 

< 

< < 

* 

< 



mm 










* 


CL 


> CO 

> 

► 

X K < 

< X 

X X 

* 

w 



CO 







o 


z 

* 

z 

3 it 

II 

X 


1— O 

4-4 X 


X 

X X 

* 



mm 

CM 







o 


4—4 

* 

1— 1 



K 3 

LL 

</> 00 

—J X 

< < 

< < 

* 

00 


f—4 

»• 







• 


o 

* 

o 

>• >- 

> 

z z 

o UJ 

-1 < 

K o o o 

* 

* 


mm 

CM 





z 


o 


3 

* 

-o 



UJ 


h- LL 

<C X o 

< 



* 

-J 


X 

*-t 





z 




u 

* 

o 

u. 


Q UJ 

h* 

4—4 

LU 

Q 



* 

< 


UJ 

\m 



o o 

z 

in 

II 

4^ 


* 


o 


Z X 

Z 

> 

o 




* 

UJ 


o 

< 

CM 


r- 

r- 

•» 

*> 



UJ 

* 

UJ 



LU |- 

UJ 

Z Q 

o z ♦ 

H- 



* 

CC 


Z 

H* 

CO 


o 

a 

4-1 

mm 

II 

z 

* 

z 

to 


CL 

X 

< Z 

LL m O 

3 



* 


mm 

4—4 

UJ 

1 


• 

• 



3 


4—1 

* 

4—4 

UJ 


UJ LU 

UJ 

< 

z 

CL 



* 

H 

*— 1 


CD 

1 

!-• 

H 

rH 

II 

II 

► 

rnm 

K 

* 

H 

H* 


O KJ 

-J 

CC 

00 00 3 

K 



* 

4-4 

•w* 

CC 


z 

1 





4—4 

f-4 

3 

* 

3 < 


z 

UJ 

UJ * 

00 1 — O 

3 



* 

O 

X 

LU 

Z 

z 


II 

II 

3 

4-4 


w 

o 

* 

O X 


•-H X 


X -J CL 4-4 LL 

O 



* 

4-4 


o 

o 

II 

II 





< 

X 

QC 

* 

a: 

4-4 


1— 1 

X > LU 





* 

-4 

-J 

UJ 

X 

II 


z z ** st 

X 

UJ 

CD 

* 

CO 

K 


UJ z 

o 

UJ O 

UJ O LU 

UJ 



* 

CL 

< h- 

X 

z 


4-4 

4—4 



X o 

3 

* 

3 

00 


X 4-4 

< 

X o 

IZQC 

X 



* 

X 

lu z a 

z z 


X 

X 

o o 

< z 

oo 

* 

00 

UJ 


H X 

UJ 

3: x 

K < < 

I— 



* 

4—1 

cc 

1-4 

o 

z z 

4-4 

oo H Q o O 

4-^ 


* 











* 














* 











* 











4* 



* 











* 














* * 

* 

* * 

* 

* * 

* 

* * 

* * * * 

* 

* 

* * * 

* 














ooooooooooooooooooooooo 


87 


in^r^ooo'0*HcvjrO'4“in>Oh-ooo>Or-ic\jfO>tin'Or-<x>0'0»Ht>jm4’u\^of s -oo 

.HrH^HrHfHCucMnjoocooocomoonjcococococorococOfOco^^Nfr^sr^'t^^ 

oooooooooooooooooooooooooooooooooo 

2Z22ZZZZZZZZZ2222ZZZZ2ZZZ22ZZZZZZZ 


00 0> O H 

sf* in in 

O o o 

ZZZ2 

^ r r ^ f ^ ^ g w w „ y _.^ 9 f ^ p ^ r t y ^ p,, ^ t pp, p,, w y ^ ^ ^ ^ ^ 

oooooooooooooooooooooooooooooooo 

333333333333333333333333 

oouoouuouooooooooouooooo 


o o o o o 

33333~53333333 

oouuoooouoooo 


O' 

o 


o 

o 


o 

o 


o 

o 


in 


X 

f— 


UJ 

o 






K 





h- 



m 




• 










+ 

















s 




— * 









X 

• 3 









O 





O 



z 













UJ 

•w 









O 





h- 







fH 









o 

X 

















II 




1 









z 










a 





o 







f-H 









fH 










o 





o 



z 




*> 









w 

> 









r- 








•H 

o 











* 

• 

a: 







in 


o 








X 

CVJ 



w 









co at 

h 









• 





z 



K 




< 









't o 








o 


fH 





f-H 



f-4 

O 



K 









• 

1 







H 


• 





X 




K 



UJ 









o ~ 










UJ 


o 



h 







CO 


O 







K O 

h- 







O 


O 


o 



e 



o 

O 



+ 


O 







• 

z 







o 


• 


rH 



UJ 

— * 


r- 

O 



-■* 


CO 







o a 

*— 1 














O 

»—» 


o 




rH 









O ul 

o 





z 


«■* 


•H 

fH 

II 



• 

•» 


• 

«■» 



1 

•H 

II 







• 

CL 





z 


H 


1 

o 



z 


00 


fH 

rH 


z 

H O 












H 


z 

M 

• 


f~4 

fH 



z 

*— i 



• 

0 


z 

•» 

00 

A 







Z 3 

U 

O 




•» 


Of 


• •* 


fH 


•* 

•«* 

< 


Of 

a 


r* 

CO 


H 







Z w 

O 



1 

a 

(\J 

II 

UJ 


00 

11 



fH 

00 

h 


UJ 

UJ 


cm 

w 

II 








Z X 


II 



fH 



• 

fH 

w* 


1 




UJ *H 


• 

• 

CM 


< 


1 


z 




rH 

• Ui 

z 


fH 

3 

+ 

II 




< 



UJ 

II 

< 

CO 

UJ 

Z 

X 


ll 

H- 

A 


Ul 

IH CO 


CO 


uj a 

o 

«■» 



3 


hh o K K 

K-# 

*-« => 


►— 

II 

3 >-« 

O h- 


UJ 

*— « 

•H 3 

X fH o 

•H 

♦ 

o z 

f-H 

3 


1 


♦H 


O CO UJ 

W 

x-" 

z 


UJ 

II 

Z X 

o «/) 

•— » 

CO 

w 

w 

z 

I- 





• rn 

K 

w 


X 


X 

-J 

UJ 

CO 

X 

X 

1— 1 


CO 

z 

f— 1 

K 

-J 

UJ 


w 

X 

X 

M 

»-H O 

II 

HQ-5T- 

o 

X 

3 Z 


in 

LU 

W 

I- 


UJ UJ 

h- h 

■o— ' 

Z M 

1- 

w 

w 

K 

O' 

w 

UJ UJ h 


h- 


h- 


w w 

Ui 

Ul 



II 


Q 




Q Q Z 



^ X 

z 






Q 

Q Z 

II 



II 

II 


-J 

o 

II 

II 


O 

Z LL 

>- u_ 

Z Z O O U_ 

X K— 

o u. 

LL 

> 

O LL 

z z o 


o 3 

o 


u u 

Ul 

z 

X 

Q 

*— i 

HH 

CO 

•H 

f— i 

•H o o 


h- h-4 

o 

IH 

M 

CO 

a 

*— 1 


O 

3 O 

HH 

3 0 3 

►-H 

Ul 

H- 1 

UJ 

LL 

X 


in o 


O' 


o 

CO 


eo co 


u o u 


88 


<m 

in 

O 

z 
»— < 
o 

-> 

o 


m>tinvOh-ooo'OHOjm>t 

ininininininino>0'0>oo 

oooooooooooo 

ZZZ 2 ZZ 2 ZZZZZ 


ooooooooooo 

ooooooooooo 

* 

o — 

Q <•"* 

• <*■* "D 

CM -d 


in^NooO'OHrgm-t^'ONQO 

vO'0'O'O'0r w Nr*“h , f N, f w f w h , h 

oooooooooooooo 

zzzzzzzzzzzzzz 

►HI— I w M W *— 4 fH fH »H *— I *— I W W H • 

oooooooooooooo 

-d-d-d-d-d-d-s-d-d-d-d-o-d-d 

oooooooooooooo 


a>OHf\jfn>tinoNQ0 

h-OOODOOCOOOCOCOGOCO 

oooooooooo 

zzzzzzzzzz 

HHHHHHHHHH 

oooooooooo 

oooooooooo 


•» 

Osl 


UJ UJ "D * 

CQ CD **X - 
sf X It 
I I ~ * * 
< O <\| 


O' H 

< < 

h- h- 


fH 

LU UJ 

a 







CD — CO 

UL 






< 

* UJ * 

* 






-* h- 

X « X 

X 






— UJ 

XXX 

X 



»H 



-> CD 

X 





in 


► w 

ft 1 

ft 



o 

fH 


in 

1 







~ + 






o 


< 

T> ~ “0 

“5 

fH 


o 

1- 


H- H 

*• “0 ► 

•» ■ 



CD 



HI Q 

O •* <M 

•H 

o 


o 


CD V* 

H O H 

fH X 

h 



CD 


ft -> 


w ^ 






— UJ 

< < < 

< CM 

o 


CM 

MM* 


-> * 

1- K K 

h- # 

CD 


X 

•M 


► cm 

UJ UJ UJ 

UJ * 



• 

fH 


CM * 

CD CD CD 

CD CD 

0 m 


H- 

X 


~ * 

V W W 

w 

z 


CD 

• 


< x 


+ 

M 


• 

H 


(\J h X 

* * * 

* 

X 


** 

-J 


CM * UJ + 


a 

oo 



- • 


* * CD ~ 

A A A 

— -> 

• 


•» 



* ~ 

H H CM 

CM •* h- 

H* 

o 

r*- 

•H 

O' 

— -> * - 

o o a 

Q vO H 

-J 

H 

w 


>fr 

-> •* X O' 

* * * 

* ~ 

* 


< 

♦ 


► H X 

XXX 

X < O 

— * 

o 

h- 


o 

O' H < 

W w — » 

~ H h- 


ft— 

UJ 

fH 

»- 

^ + H- 

sss 

UJ 

H 


CD 

w 

O' 

< < UJ CM 

CD CD CD 

CD CD O 

1 

o 

*•»* 

X 

<4- O 

hi CO O 


CD 

M 

CD 

• 

w 

o 


I + + ft + 


-5 T") -5 n H 


>o 

H 

o 


<w- 

< 


o 

r- 

Q 


oc 

o 


cc 

o 


o ** 
cd x 


M M 

h- O * 


h i\i 4 in ro 

oo 


H 



UJ 



• 


X 


CM 







00 

T> -5 

UJ • * 


w» w wo 

w 

w 

• 





CD 



H 


• 


X 






A 

• 


CO CM X 


< < < < 

< 

UJ M 


O 



+ 



• LL 


H 


• 




CM 


o 

h- 

O CM 

1 X 


h- h~ h- 

h~ 

1— 

CD -D 


CD 



M 

o 


O UJ V 


-J 


h 




•H 


• 

CD fH 

•H fH 

1 — -f 


UJ UJ UJ 

UJ 

UJ 

• f 




Z fH 

o 


UJSi- 


• 

z 

CD 




•> 


UJ 

• *• 

W 

-5 M 


CD QQ CD CD 

CO 

m m 


M 


z 

1 

|H CO 


• M -0 


0 m 

z 

• 

o 



fH 

z 

z 

in 

< < 

M ► -5 







fH 


z 

M 

o 


*m m 


“D Z 

a 

o 


H 

fH 

z 

• 

fH W 

►— h- 

“DO* 


II tl II 

II 

II 

* < 


• 


►> 

•» 

CO II 


“> * fH 


•» 

•* 

•H M- 



f O 

#• 

0 m 

* < 

UJ UJ 

► H H 





in x 


a 


CM CO 



* 0* fH 


h- fH 




+ 

fH f*— 

fH 

fH 

in x 

CD CD H w *H 


M M -M 

0 m 

M 

^ X 


UJ 



w 

II — • 


00 w w 


w 


4 

11 



fH Q 



w JT 


< w 

X 

J 

J 

J 

-D 

“D 

< < 

*D 

• CO 

II 

< 

fH 


w < < < 

II 




X 

• 

II 

X 

< < 

* * 

<h < 


•* •> f 

f* 

•» 

X CD 


X 



H* 

— ft 

UJ 

< h h 


I- 


fH 


UJ 

“D 

fH 

UJ 

X CD 


hlilh 

\HMfO 

'fr 

in 

X 

It 

O h- 

•H 

UJ 

fH fH 

DhlUUJhUJ 

fH 

w 

fH 3 

•H 

X 

fH 

o 

X 

UJ UL UJ CD UJ 


WWW 


w 

< II 


O v/> 


CD 

w w 

Z UJ CD CD to CD 


X 

w 

z 


“0 it 


z 

< M 


CO ^ CD CD < < < 

< 

< 

CD 

z 

J UJ 4 

w 

X X 

fH 

CD 

UJ 

w 

in 

m* 

X 

fH 

II 

•H 

O' fH 

CD 

II II 

•8' 


XXX 

X 

X 

— Z 

*H 

w 

K 

fH 

w 

UJ UJ K 

— II II 

H 


H 

w 

UJ H 


-Z4 

'W 

— Z 


II II X 

II 

XXX 

X 

X 

HH X 





O O Z 






O Z X 

«H 



M 

fH CM 

X 


< <t < 

< 

<r 

U X 00 

u. 

> O LL 

ZZOU.r-tCMVU.OU. 

Z O -5 U. X o u. 

LL X 

O O CD X 

X 

CD CD CD 

CD 

CD 

♦H l/) 

fH 

fH 

CO o 

►H 

fH fH O 

•H X X 

CD 

HI 

O IH> 

•H 

O «H 

•H CO Q fH 

•H 00 


in o oo 

•h fH *4" 


o 


o 



89 


0'Q*M<Mc^'4"U^>or^ooa>OrHc\jmNtinvGh-cQ 

ooa'oc'o^oa^o^a'cho'ooooooooo 

OOOOOOOOOOOr-tr-4rHrH*"lr-l*Hf««rH 

2Z2Z2222222222222 2 22 

hhhhhnhhhnhhhhhhhhhn 

aooooooooooooooooooo 

oooooouuooouoooooouu 


<\J 

in 


in 






CO 






O 






O 

in 






h 


2 



2 

H 









2 



2 


o 






O 


2 




O 

H 






CD 





II 

o 

O 

2 

2 

2 

2 

2 



II 



2 


CD 

2 

2 

2 

2 

2 

2 


2 



1—4 

2 


2 

2 

2 

2 

2 

1-4 


4—4 



2 

4-» 

A 

Ok 

•» 

*■ 

*> 

*• 

X 


X 



h— 

2 

o 

H 

(M 


in 

00 

00 


00 



4-H 

h 

o 


w 

W 

— *- 


t 


4-4 




• 

1-4 

< 

< 

< 

< 

< 

1- -* 





—* 

LU 

• 


h* 


H 

1- 

CD 2 


— * 



o 

O 

LU 

LU 

UJ 

UJ 

UJ 

LU 

• 2 


o 



r- 

« 

CD 

00 

00 

00 

GO 

CO 

— 2 

o 

r- 



o 

•-* 

• 






2 

o 

o 



• 

2 

a 

II 

II 

tl 

II 

II o 

2 in 

in 

• 



vH 

2 

2 






2 w 


»—i 



• 

2 

2 

<•— » 

,«% 

—V. 

«— 

-* II 

- < 

ii 

• 



a 

•» 

2 

2 

2 

2 

2 

2 

m x z 


a 



LU 

00 


2 

2 

2 

2 

2 ^ 

~ X 2 


UJ 

4— • 


• 

w 

X 

2 

2 

2 

2 

2 2 < < 2 

in 2 

• 


LLi 

2 

< 

LU 

•* 


•> 


*■ z x o 

2 

2 

II 

D 

4—4 

»- 

o 

i-4 

<\J 

m 


>A Z X II 

H- 2: 

4—4 


2 2 

UJ 

2 

kw 

w 

w 

w 

*— > W 

< II 

2 00 — 

X 2 

2 

•— i 

1- 

00 

►— « 

< 

< 

< 

< 

< X O 2 CL LU X 

00 CL 

4—4 

K 

w 

•<— > 

w 

2 

X 

X 

X 

X UJ 

Z H o h UJ 

w 3 

2 

2 




2 

2 

X 

X 

x a 

►-H X 

I— o 

H 

00 o 

LL 

LL 

LL 

< 

< 

< 

c 

<2LLX00LU>2 

LL UJ 

4—1 

O 

1—4 

i— i 

CD 

CD 

o 

o 

CD 4-4 

H l/) M 

(L CO m 

4-4 CL 


o 

z 

LU 


O' 


o 

m 


in 


c\j m 
in in 


o 



SUBROUTINE OUTPUT PRINTS THE COMPLETE STATISTICAL ANALYSIS OF SJOIN AND 
C JOIN OF A SUBSET T AND D» EACH OF LENGTH L. INDEX IS THE VECTOR CONTAIN- 


90 


* * * * 


* 

* 

* 

o u o o u 


* 

* 

* 

* 

cjrosfm.or^coa' 

O rH CM ro >4“ in 

QOO>OH(M^ 4- in >ON 




* 

OOOOOOOOr-l.-ltH.-l.-CrH 

MHHH(MCM(M(MCM(MCMCM 

1 



* 

oooooooooooooo 

oooooooooooo 

X 



* 

(-(-•-H-KHH-K 

fr- 

H H h 

fr- 

1- 

fr- fr— 

h- 

fr- fr- 

fr- fr— fr- fr— fr— fr- h- 

3 



* 

oooooooooooooo 


z 



* 

v 

oooooooooooo 

o o 

oooooooooooo 

UJ 



* 










X 



* 








CO 


fr- 



* 








DC 





* 








• 


iz) 



* 








•» 


mi 



* 








X 





* 








H 





* 



< 





H 


LL 



* 



in 





•* 


M 



* 



•* 





w 


X 



* 



X 





CM 


CO 

• 





o 





CO 


•— i 

fr- 


* 



f-4 





m 



ill 


* 



•* 





*• 



CO 


* 



CM 





X 


• 

CO 


* 



Mi 





CM 


CC 

3 


* 



•» 





H 


UJ 

CO 


* 



m 





*• 


GO 



* 








m 


X 

UJ 


* 



1 





(M 

CM 

3 

X 


* 








< 

Z 

Z 

fr- 


* 



m 





m 

» 




* 



•* 





#* 

rH 

Z 

o 


* 



CM 





X 

z 

Z> 

z 


* 



w 





Ml 

*> 

cC 

1— » 


* 



•» 





Ml 

A 


Q 


* 



m 



_J ' 


•» 

00 

UJ 

UJ 


* 

A 





** . 


m 

•» 

X 

UJ 


* 

CM 


-J 



f— 1 


9-4 

H 

fr- 

o 


* 

CM v> 


< 





to 



UJ 


* 

•* co 


> 



II 


CO 

II 

CO 

C* 


* 

in v. 


CC 





a: 


Mi 

a 


* 

~ h- 


UJ 



Mi 


m 

Mi 




* 

< z 


fr- 



** 


•> 

•» 

3 

CO 


* 

X m. 


z 





X 


M 

fr— 


* 

X CC 


Mi 



Ml 


Ml 

3 


z 


* 

< CL 


■» 



W M 


Ml 


• 

UJ 


* 

O Z — 


** 



h 


•» 

Mi 

Z 

X 



•» •* ^ 


X 



H X 


• 

o ^ 

M 

UJ 


* 

— ^ <M 


O CM 



— CM 


H 

H < 

o 

-J 


* 

cn cm w 


m* m — * 



M« M 


GO 

fr— 

-o 

UJ 


* 

CM ^ ^ 


•'CM- 



w ro 


• 

O UJ 

u 



* 

► X Mi 


(Mm* 



o • 


•» 

O CO 


o 


* 

00 o * 


M M 



» |M 


X 

• . 

X 



# 

w Z — , 


•* m 



— LL 


CM — 

*. 

o 

Q 


* 

0 3 ^ 


9-4 UJ 


fr- 

M4 M 


Ml • 

O O CM 

Qt 

z 


* 

fr- Q- CM < 


Hi X 


LL 

w M 


H fr- 

r- z 

LL 

< 


* 

to z X 


♦* a: w 


Ml 

X • 


•' m 

II o •* 




* 

•* •» •—'X 


-)UJh 


X 

Ml CO 


9-4 H 

fr- fr— • 9-4 

00 

h- 


* 

— V. H< w 


H CO 


CO 

W LL 


< X 

LL LL CM Ml Z 

UJ 



* 

ro «— i •» O h* 


X 1 


Mi 

► CM 

- CO 

Mi Mi Z • 

O 

o 


* 

CM V* -I * •* •* 


— 3 



— CM O 

•* H 

XX ar-M 

O 

U1 


* 

•* O — 1 * 9-4 O O fr— 

fr— 

CM Z UJ 


+ 

(O Mi Ml X H 

to co m. uj in 

U 

a 


* 

(M<in^w|-WUL 

u. 

► o 



** *• 

•* 

^ • 

Mi Mi M* • *» 


QC 


* 

HUJ-HXV)\h 

HH H Z Z -J 


fr— • 

fr- r-4 z 

Ml — H* 

K 

< 


* 

~ a: a ~ uj * i 

X 

Z 3 < 

*• 

Ml 

z 

z 

•* Mi 

X 4- *f 1 -0 Z 

00 

o 


* 

<zwoo<ujcoco 

H-* CC fr- 

H 

w 

Mi • 

Mi 

• o 

f x •* Mi 

UJ 

CO 


* 

fr- ► Z h flO hi 

Mi 

OC O CO 


MM 

CL W 

oc o -> 

Ml Ml CM w oo OG 

K 

Ml 


* 

UJ «4“ « M 111 < 


CL » Mi 

II 

►M 

CL O CL 

• m 

Ml • M< * 0, 


a 


* 

CO * * m» CQ J + 

+ 

Z w Q 



z ~ 

z 

w *> 

II -fr -fr UJ < z 

UJ 



* 

(* OG w N* 


w 

NH 

II 

w 


X 

1 O fr— — 

X 

LL 


* 

COUJOJhUJZZH 

-J 

fr— 



fr- 


fr— r4 

-0 3-0 t LU 

fr- 

o 


* 

* o o o a o 


UJ < 

o 

«■» 

UJ < 

UJ < <H J ”3 CD UJ 




* 

J UJ UJ J Ul X X II 

II 

fr- x m 

Ml 

Ml 

I- X 

fr- 

X H 

4* « II w W fr- 

o 

CC 


# 

< fr- fr- < H X X 


M OC w 



h o£ m ac • 

II »-< 

z 

Ul 


* 

UJZZtUZOOHNdCOOO^ 

of OaCON 

O Ml CM LL IL CC 

*•4 

GO 


* 

Q^HWOCwUO^ 

1— 1 

3: LL 

o 

Mi 

111 JL. QZZwm J 




* 



X 





X 





* 



CM 


o 

CO 


CM 





4fr 





9-4 



O 


* 

* 

tf 

* 








Ml 


o 

c 

c 

c 

c 











91 


oooOi-i<>Jcr»^mor^oochOr-*<Mm<ftnor^coo*o 

ooooooooooooooooooooooo 

KHI-f-H-Kf-KI-f-HKKKH-KH-H-l-KI-K*- 

O O o O 0.0 ooooooooooooooooo 


X 

00 


GO 


o 



_J 



H 




U1 



• 




> 



to 




• 



• 




H 



H 


“5 


X 



X 




o> 



CO 


X 


H 



r-4 


UJ 


m 



H 


a 


OC • 


CM 

» 


2 


CC — > 


2 

CM 


r— i 


O CM 


H 

GO 


•* 


o co 


H 

m 


CO 


• w 


2 

H 


2 


* O 


H 

X 


H 


X _J 


0 m 

CO 




O' UJ 


00 

r-4 


in 


H > 


H 

H 


H 


• • 


H 

• 


r-4 h 


M H 



CM 


O 


H X 


II 

< 


II r-4 


co <r 



• 


H4 


<W H 


1— « 

H 


H H 

M. 

o • 


H 

X 


h CM 

CM 

H* a: 


A 

CO 


CM -* H4 

HI 

00 0 C 


r- *> 

r-4 


r^ ► 

CM 

- o 


H 

H 


H X 

H 

► o 


o ~ 

m 


O H4 r-4 

X 

X » 



H 


h- •* 

•-4 

CO H 


Q 

CD 


< in 

*» 

H X 


O f~ 

m 


o x • 

in 

• O' 


O 00 

H 


O x in 

• 

H 


— * . 

X 


C hH 

in 

r-4 • 


— * H 

CO 


— O Cl 

*-4 

< — 


O O CM 

H 


o in 

o 

^ CM 


r- 2 

H 


O - H 

00 

O 03 


no- 

H 


• CO X 

H 

1- ^ 

h H 

• f-i 

r-4 


K O 2 r-4 

X 

— 00 O 

a. u 

CMHZ 

— < 


U • H 

H 

CO • 

HI •—« 

2 • 

- 


m a-cg 

H 

O H 00 

X X 

Of ^ 

O - 


X UJ CO H* 

CM 

H X • 

00 to 

uj in 

r-4 X 


to • r-4 H 

HI 

h H •* 

1-4 H-4 

**% • H 

* CO 

X 

M — *« X 

CM 

K r-4 X 


r-4 — H- 

h H 

X 

-5 h CM 

H 

2 h CO X 

+ + 

1 -5 2 

2 <► 

H 

4- -2 r 

• 

M * HO* 


X ► H 

H4 » 

r-4 

in m * 


a: o • h h m 

wOOOC 

OC O CM 


r-4 — a: 

m 

G- • • 


• w a 

a • 

II 

<X Q. • 

^OJZ^N II 


UJ < 2 UJ 

2 — 1 


4 XZ^UJ 

o 

— < 




*•> 

X ^ z> 

H 2 

K ^ -0 

-> ■> 

• UJ 2 

~ h- -J 


-> < 1—2 2 

<t M 

UJ < O 


-)C0 111 M 

UJ < 

CM 

O u» < #-h a: 

x »- k x h- r- 

II II 

W H K H X II 

r-4 

II K X H 3 

CC 2 

M (/) 


*-« 2 

a: 


wflC ZhQ 

OOotO - 0«-«<M 

uuu.a:ootO^Orou.a:oouj2 

LL O 3C LL O 

2 2 

M M i u J IL X 

OZHlltUOnil 


X 






in 4- 

CO 





fO CM 


o 



O 


f-( r-l 


H 



r-4 





SUBROUTINE ABYSS (N,NPTS,S,DELTA f ZETAI ABYSSOOl 


92 


* 

* 

* 

* 

* 

* 

* 

* * 

* 


* 

* * 

* 

* 

* 

cm co ^ in 

o 

b- 

CD 

O'OHcgm 

>t in 

no r- 

CO 

* 







00 







* 

O O o o o o 


fH 

rH 

rH 

* 







LL | 


1 





* 

o o o o o 

o o 

oooooooooo 

* 


> 





uj Z 


z 





* 

to to oo 

OO 

to 

to 

CO 

to to oo 

(S) 00 

00 

00 

00 oo 

OO 

* 


< 





o o 


o 


LU 

LU 



to to oo to 

O0 

00 

to 

to 00 to 

O0 00 

00 00 

to oo 

00 

* 


cl 





o o 


o 


00 O 

o 


* 

> >• >- 

> 

> 

>• 

> 

> V > 

> > 

> 

> 

> 

>- 

> 

* 


cl 









LL Z 

z 


* 

CD CO CO 

CD 

CD 

CD 

CO 

CD CD CD 

to CD 

CD GO 

CD 

CO 

CD 

* 


<t 





Z cL 


a: 


UJ UJ 

LU 


* 

< < < 

< 

< 

< 

< 

< <■ < 

< < < 

< 

< < < 

* 







C UJ 


UJ 


o o 

o 


* 













* 

' 

LU 





LU X 


CL 


u •-* 



* 













* 


X 





X O 


CL 


UL 

LL 


* 













* 


I- 





wl 


3 


z z 

Z 


* 













* 







X 




< o 

o 


* 













* 


X 





1- UJ 


UJ 


UJ U O 


* 













* 


H- 





1 X 


X 


X 



* 













# 


1-4 





HH ►— 


h 


QC 

CL 


* 













* 


X 









X UJ 

UJ 


* 













* 







Z Z 


Z 


H X 

CL 


* 













* 


to 





o o 


o 


• o 

CL 


* 













* 


LU 









i-i _j 

3 


* 













* 


O 





Q O 


Q 





* 













* 


< 





UJ UJ 


LU 


z z 

Z 


* 













# 


U. 





oo to 


OO 


o o o 


* 













* 


QC 





< <t 


< 





* 













* 


LU 





GO CD 


CD 


a o 

a 


* 













* 


h— 









UJ UJ 

LU 


* 













♦ 


Z 





DC oC 


CL 


to 00 

to 


* 













* 


i— « 



mm 


LU UJ 


LU 


< < 

< 


* 













* 





z 


> > 


> 


CD CD 

CD 


* 













* 


UJ 



•* 


< <f 


< 





* 













* 


X 



vO 


-1 -J 


-J 


CL CL 

0C 


* 













* 


K 



w 






UJ UJ 

LU 


* 













* 





< 


X X 


X 


>» > 

> 


* 













* 


o 





H* H- 


H 


< c 

< 


* 













* 


b- 



-J 


1 » 


1 


-J -J 

~J 


* 













# 





LU 


»— 1 4—4 


*-4 





* 













* 


00 



o 



K 


K 

X X 

X 


# 

CL 












* 


X 





U. UL 

M 

LL 

M 

*- b~ 



* 

or 












* 


H 



» 


o a z 

O 

X 

1 • 

1 



00 












* 


CL 



1 



*-4 


»— 4 

l-H *— » 

4-H 


* 

o 












* 


Ul 




cL 

oo oo 

-J 

00 

— J 




# 

•» 












* 


o 



00 O 

to to 


to 


U- LL 

LL 


* 

mm 












* 





cL 

b- 

UJ UJ 

UJ 

UJ 

UJ 

O O 

o 


* 

>fr 












* 


LU 



LU 

u 

z z o z 

o 




* 

CM 












# 


X 



> 

UJ 

2* ^ 

z 


Z 

X X 

X 


* 













* 


H- 



< 

> 

O O LU 

o 

UJ 

h- h 

f— 


* 

to 












* 





-J 


H-4 H-4 

a 

4—1 

a 

a. a 

CL 


* 

•* 












* 


00 




K 

X X 


X 


UJ UJ 

UJ 


* 

■mm 












* 


UJ 



UL 

1 

1- K 

U- 

H LL 

a o 

O 


* 

H 












* 


h- 



O CO 








* 

<-» •« 












* 


< 




1- 

II II 


II 


ii ii 

II 


* 

M sO 












* 


-J 



CL 

z 








* 

1 — 












* 


3 

0 


UJ 

UJ 

«nt mm 


mm 


mm mm 

•*—. 


* 

o < 












* 


O 

00 


CD O 

»— 1 1-4 


4—4 


1-4 1— I 

1-4 


* 

•“ h" 












* 


-i 

h~ 


X 

3 

•> •» 


•* 


► •* 

•» 


* 

X UJ 








CM 

CM 


X 

* 


< 

Oo 




CO 


sfr in 

o 


* 

1 KJ ^ 











X 

* 


u 

Z 


z 

00 

w w 




w w 

w 


* 

<t •‘CO 


Z 






1 

O 


X 

# 







< < 


< 


<t < 

< 


* 

W — X. 


oc 







H 



* 


00 

a: 


M 

II 

b~ 


H- 


h- f— 

1- 


* 

rH h- 


3 



O 



mm 



— 

* 


oo O 




UJ UJ 


UJ 


UJ LU 

LU 


* 

00 ^ Z 

a 

H 



O 



rH 

O 


rH 

* 


> 



Z oo 

M hsl 


M 


ivj M 

M 


* 

*>Ow 

fH 

UJ 



z • 

z 



O 


»• 

* 


CD U 











* 

-J w cc 

«w 




Z ' 0 o z 


+ 



1- 

* 


< 

UJ 











* 

< < CL 

00 




•» •* 

m 





z 

* 



> 











* 

uj t- Z ►- 

a 



HH II 

w4 


X 

rH 

•H 

4—4 

* 


UJ 












* 

a: -j 

a 

i-i 





vH 

w 

• 


oc 

* 


Z 

a 











* 

uj sf Z 

• 



II II ~ 

II 


00 

UJ 

4* 

CL 

* 


*— 4 

z 











* 

H- O * 


a 

• 

o 

•—4 


1 

H 

O 


Z 

* 


H- < 











* 

»-h CL 

cL 

LU 


Q 

I 

J 

t 

X 


Q. 

• 

X 


* 


3 












* 

U oo UJ 

UJ 

• 

z 

• 

3 


X 

Z 

3 



* 


O 

< 











* 

w# OOZ 


o 

in in w 

o o 



3 

II 

UJ 

* 


a: 

H- 











* 

J J LU 

uj 

w 

II 


CM CM < 

CO 

II 

II 

w 


h- 

* 


CD 

-J 











* 

Q. < h 

b- 



u 

K 

II 




X 

1-4 

* 


3 UJ 











* 

Z UJ Z Z LL 

z 


o a id 

o 

X 

3 

UL 

X 

CL 

* 


00 O 











* 

»“4 QL *-» 

•-H 

1— • 

ZOOQN 

a 3 

X 

3 

4—1 

X 

X 

* 














* 













* 














* 






in 







# 














* 






CM 







* 

* 

* 

* 

* 

* 

* 

* * 

* 

* 

* 

* * 

* 

* 

* 














ooooooooooooooooooo 


93 


o> o »-• cm co ^in>ohQOOOH(Mm>j , iA>o>ooo'OH(Mro>tiA'Of^ 
p-*fM<McM<M<MCMCMfMc\jc\icocorncocomcOrOcoco>4“M’<4'frM ,, ^>t4’ 
oooooooooooooooooooo o o o o o o o o o 

COl/)CO</)CQCOCOCOt/)COlOCOCOcOCOoOl/)i/></)COCOt/)COCOCOcOCOLOCO 

tOOOOOtOU)(^OOlOOOU5lOl^OOt0^tOUOt/)^OOU)OOOOOO^OO^CO(/) 

COGOCOcOcOCOaOGQCDCOCOcQCQcOCOCOCDcOcOCOCOcOCQCOCOaQCOCOcO 

<<<<<<<<<<<<<<<•!<•<<<<<<<<<<<< 


CO 









4- 









Z 









*-* 









o 









a 


sD 







< 


O 

o 






4- 


4- 




Mft 



< 



o 


rC 

r-4 

H H 


O 


o 

H 








o 



4 

4 

4 4 


CO 



O 









o 


X 

X 

X X 


z 





#• 

p» 

•» •* A 

M> 

< 





CM 4* 

CM4I 

X 

X 


ME* 



w 


W W »> 

•» 

1- 


o 

o 


< 

< 

< < CM 

CO 



Q 

o 


H 

h 

4- 4- ~ 

w 

a: 


0 

• 


-J -J 

J J< 

< 

LU 


O -> 

o — 


IU UJ 

UJUJh 

H 

X 


• 'CM 

• CM 


o a 

O O LU 

UJ 

UJ 


LU * 

UJ * 




NJ 

isl 

a. 


* 

-» * 


♦ 

# 

♦ * 




« CO 

• CD 




4 

4 

</> 


«*» 

’ 



Mi 

0 m «m 


<t 


CM 1 

<M f 



-> 

— 

a 

X 


* 

♦ 

QC 


-5 

-> '■o X 

X 



* CM 

♦ CM 


w 

w 

w W •» 

»• 



00 * ~ 

CD * 


CO 

CO 

CO CO rH 

fM 



# QC 

♦ 

~ 4- 





<M 


• — * 

1 M 

1 

4* 

4 

1 • < 

< 

1-4 


X O 

t— » 

► O O 



H 

4- 

_ •* 

A a 

CM •* O 

CM 

HOIO 

Mi 


— ui 

LU 

• 

H H 

♦ mo* 

♦ ro 


•-4 

r4 — 

f-i t* CM 

Nl 


QC 

4 4 


* ~ H 
-> < 

CM 

4 ^ < (\l 
— < t- 

o X 

(M + + *• 

4 4 1 | 

UJ 

> 

X X 


X K O N* 
_J O <L 

M|-UJ< 
X •* J Kl w 

► H 

O X X ^ 

o x x o a 

< 

-J 

•i •» 

^ CO 


CO UJ 

^ o 

II 

X CO UJ 
* ~ 0 * 

u 

H fO h 

H rH CO II II 

• 

W ^ < w (M 

<<p4Hh • 


-4 < w 

H 4- Q- 


^Vvlil 
CM < < Nl 

W‘ w 

'*■* < ^ *■* 


H H 


JQCUJ 

r-* O 

M -J a: 

X 

H •— 

X K 4- X X 


-J -J 

4 

uiao 

•» o o 

LU o 4 

•* 

1 -1 -» 4 

•CO J J * CO •* 

h H UJ UJ 


O CO • 

-5 • 

w o CO 

-o 

UJ UJ 

^ uj UJ in oz 

< 

O O 

-o 

W o X 

W o o 

-Oh 

w~ 

-5000 

wOOO — Owcx 

X 

1) 


mt V 

<f h~ 

m ~ 

< 

w 

< 4- < h O 

CC 

It It 

11 

It 

K It 

n it 

4- 

It It II 

t- II II 4- 4-4-0 

o 



LL a. 

UJ O 

O LL 

LU 

LL 

UJ O UJ O UJ LU z 

LL 

"7) < CO 

-7 

QC *-* 

l 

G 

T 

D 

I 

p 

T 

M 

H<flOQKio<coNONa:uJ 

i-4 

CM 

CO 


>3* 

UN 


o o 

o o 


-» CM m 


94 


APPENDIX B. TABULATION OF PREFERRED REGRESSION FITS 


Table 6. Tabulation of preferred regression fits 
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7-04E 349.8 1717(11) 2396(12) 67.4 90.4 47.8 
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Table 6. (Continued) 





10- 3W 370.1 1879(11) 2267(08) 76.5 215.2 13.0 

11- 1E 357.7 1936(06) 2133(02) 14.4 - - 5312(14) 94.5 99.7 89.4 

11-1E 356.6 2003(07) 5433(15) 89.7 95.3 84.2 

11-1W 367.7 1804(06) 4776(15) 68.6 73.1 64.2 

11-1W 367.7 1753(12) 4965(10) 73.9 75.9 72.0 
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